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^ ■ ABSTRACT 

Context. Towards the end of the evolutionary stage of the Asymptotic Giant Branch (AGB) the atmospheres of evolved red giants are 
["T | , considerably influenced by radial pulsations of the stellar interiors and developing stellar winds. The resulting complex velocity fields 

severely affect molecular line profiles (shapes, time-dependent shifts in wavelength, multiple components) observable in near-infrared 
spectra of long period variables. Time-series high-resolution spectroscopy allows us to probe the atmospheric kinematics and thereby 
study the mass loss process. 

Aims. With the help of model calculations the complex line formation process in AGB atmospheres was explored with the focus 
on velocity effects. Furthermore, we aimed for atmospheric models which are able to quantitatively reproduce line profile variations 
found in observed spectra of pulsating late-type giants. 
Methods. Models describing pulsation-enhanced dust-driven winds were used to compute synthetic spectra under the assumptions 
of chemical equilibrium and LTE. For this purpose, we used molecular data from line lists for the considered species and solved the 
radiative transfer in spherical geometry including the effects of velocity fields. Radial velocities (RV) derived from Doppler-shifted 
(components of) synthetic line profiles provide information on the gas velocities in the line-forming region of the spectral features. In 
addition, we made use of radial optical depth distributions to give estimates for the layers where lines are formed and to illustrate the 
effects of velocities in the line formation process. 

Results. Assuming uniform gas velocities for all depth points of an atmospheric model we estimated the conversion factor between 
gas velocities and measured RVs to p = u gas / RV ~ 1 .2-1.5. On the basis of dynamic model atmospheres and by applying our spectral 
synthesis codes we investigated in detail the finding that various molecular features in AGB spectra originate at different geometrical 
depths of the very extended atmospheres of these stars. We show that the models are able to quantitatively reproduce the characteristic 
line profile variations of lines sampling the deep photosphere (CO Av = 3, CN) of Mira variables and the corresponding discontinuous, 
S-shaped RV curve. The global velocity fields (traced by different features) of typical long-period variables are also realistically 
reproduced. Possible reasons for discrepancies concerning other modelling results (e.g. CO Av = 2 lines) are outlined. In addition, we 
present a model showing variations of CO Av = 3 line profiles comparable to observed spectra of semiregular variables and discuss 
that the non-occurence of line doubling in these objects may be due to a density effect. 
Conclusions. The results of our line profile modelling are another indication that the dynamic models studied here are approaching a 
realistic representation of the outer layers of AGB stars with or without mass loss. 
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1 . Introduction visual and infrared wavelengths (IR; e.g. Lancon & Wood l2000l 

Gautschy-Loidl et al. l2004l GH04, Aringer et al. 120091 . 



Stars on the Asymptotic Giant Branch (AGB) represent objects _ , . , , __, , , . 

£1 . . . f . . . f n o o it \ • On the upper part or the AGB, the stars become instable to 

or low to intermediate main sequence mass (xO.ti-tS M B ) in a ,. , r± \ *\ _ . , ' 

, . , .. , Yin.-! i , a .■ . strong radial pulsations. This leads to a pronounced variability or 

late evolutionary phase. While they exhibit low effective tem- , . , „ • > >• r , • , • 

/^cn/w\ .l • i • v u i e the emitted flux with amplitudes or up to several magnitudes in 

peratures (<3500K), their luminosities can reach values of up to . . ,, . . K, „. , ~„„^ „. , . . 

c m4 t c *t. a/->t> 1 ■ the visual (e.g. Lattanzio & Wood 2004 ). Since the variations oc- 

a few 10 L Q at the tip of the AGB, placing them in the upper , • , r r , , ,™ , 

• , . c xx . D r, j . cur on long time scales of a few 10 to several 100 days, pulsating 

right corner of the Hertzsprung-Russell diagram. Compared to & . . , , . , . ,r 

.. . ■ e . . e . . i e AGB stars are often referred to as lone period variables (LPVs). 

the atmospheres of most other types of stars the outer layers of . • • , , 

, j j . . , i ui i a. In the past, different types of LPVs were empirically classified 

these evolved red giants have remarkable properties (see e.g. the , , . . , ,. , , r / . 

i f-\ . c e Tj -e fwvti u l r^un/i\ according to the regularity of the light change and the visual 

review by Gustafsson & Hofner 2004 below GH04). ,. , . , at, ™ 

light amplitude: Mira variables (regular, Ak>2.5 ), semiregular 

In the cool and very extended atmospheres (extensions of variables (SRVs, poor regularity, AV<2.5 m ), and irregular vari- 

the same order as the radii of the stars; up to a few 100 R ), ables (irregular, AV<l-2 m ). Major advances in our understand- 

molecules can form. Their large number of internal degrees of ing of the pulsation of AGB stars were achieved by exploiting the 

freedom results in a plethora of spectral lines. Thus, molecules data sets of surveys for microlensing events (MACHO, OGLE, 

significantly affect the spectral appearance of late-type giants at EROS), which produced a substantial number of high-quality 



2 



W. Nowotny et al.: Line formation in AGB atmospheres including velocity effects 



lightcurves for red variables as a by-product. According to the 
pioneering work in this field by Wood et al. ( 1999) and Wood 
(2000), and to a number of subsequent studies (e.g. Lebzelter et 
al. I2002cl Ita et al. I2004al I2004bb it is probably more adequate 
to characterise LPVs according to their pulsation mode than to 
their light change in the visual as it was done historically. From 
observational studies during recent years (see e.g. Lattanzio & 
Wood 120041 ) it appears that stars start to pulsate (as SRVs) in 
the second/third overtone mode (corresponding to sequence A in 
Fig. 1 of Wood 2000) and switch then to the first overtone mode 
(sequence B). Light amplitudes are increasing while the stars 
evolve and finally become Miras. In this stadium they pulsate 
in the fundamental mode (sequence C) and show highly periodic 
light changes. Observational evidence for this evolution scenario 
was found for LPVs in the globular cluster 47 Tuc by Lebzelter 
et al. d2005bl ) and Lebzelter & Wood (I2005T ). 



The pulsating stellar interior of an AGB star severely influ- 
ences the outer layers. The atmospheric structure is periodically 
modulated, and in the wake of the emerging shock waves dust 
condensation can take place. Radiation pressure on the newly 
formed dust grains (at least in the C-rich case, cf. Sect. l2.U leads 
to the development of a rather slow (terminal velocities of max. 
30km s~') but dense stellar wind with high mass loss rates (from 
a few 1(T 8 up to 10~ 4 M yr" 1 ; e.g. Olofsson 12004] ) . 



As a consequence of these dynamic processes - pulsation 
and mass loss - the atmospheres of evolved AGB stars eventu- 
ally become even more extended than non-pulsating red giants in 
earlier evolutionary stages. The resulting atmospheric structure 
strongly deviates from a hydrostatic configuration and shows 
temporal variations on global and local scales (Sect.[2]i. The 
complex, non-monotonic velocity fields with relative macro- 
scopic motions of the order of lOkms -1 have substantial in- 
fluence on the shapes of individual spectral lines (Doppler ef- 
fect). Observational studies have demonstrated that time series 
high-resolution spectroscopy in the near IR (where AGB stars 
are bright and well observable) is a valuable tool to study at- 
mospheric kinematics throughout the outer layers of pulsating 
and mass-losing red giants (e.g. Hinkle et al. 1 19821 from now on 
HHR82, or Alvarez et al. 120001 . Radial velocities (RV) derived 
from Doppler-shifts of various spectral lines provide clues on 
the gas velocities in the line-forming regions of the respective 
features. A detailed review on studies of line profile variations 
for AGB stars (observations and modelling) can be found for 
example in Nowotny (2005 below N05). 



In two previous papers (Nowotny et al. 2005 a + 2005b 
from now on Paper I and Paper II, respectively) we investigated 
whether observed variations of line profiles can be compre- 
hended with state-of-the-art dynamic model atmospheres. We 
were able to show that the used models allow to qualitatively 
reproduce the behaviour of spectral lines originating in differ- 
ent regions of the extended atmospheres. The work presented 
here can be regarded as an extension of the previous two papers 
about line profile modelling. The aim is to shed light on the in- 
tricate line formation process within the atmospheres of evolved 
red giants with an emphasis on the velocity effects (Sects. [3]|4]|5] 
and 17. 41 . In addition, we report on our efforts to achieve realistic 
models, which are able to reproduce line profile variations and 
the derived RVs even quantitatively (Sects. |6]and|7]). 



Table 1. Characteristics of the dynamic atmospheric models for 
pulsating, C-rich AGB stars (for a detailed description see text 
and DMA3) used for the modelling. 



Model: 


W 


S 


M 


L* 


[L e ] 


7000 


10000 


7000 


Mi, 


[M e ] 


1.0 


1.0 


1.5 


T* 


[K] 


2800 


2600 


2600 


[Fe/H] 


[dex] 


0.0 


0.0 


0.0 


C/O 


by number 


1.4 


1.4 


1.4 


R* 


[Re] 


355 


493 


412 




[AU] 


1.65 


2.29 


1.92 


log 8* 




-0.66 


-0.94 


-0.61 


P 


[d] 


390 


490 


490 


At( p 


[kms- 1 ] 


2 


4 


6 


k 




1.0 


2.0 


1.5 


Am bo | 


[mag] 


0.21 


0.86 


1.07 


(M) 


[Moyr 1 ] 




4.3- 10- 6 


2.5-10- 6 


(u) 


[kms" 1 ] 




15 


7.5 


(fc) 






0.28 


0.40 



Notes: Listed are (i) parameters of the hydrostatic initial model, (ii) 
quantities derivable from these parameters, (iii) attributes of the inner 
boundary (piston) used to simulate the pulsating stellar interior as well 
as the resulting bolometric amplitude Am bol , and (iv) properties of the 
resulting wind. The notation was adopted from previous papers (DMA3, 
DMA4): P, Au p - period and velocity amplitude of the piston at the in- 
ner boundary; f L - free parameter to adjust the luminosity amplitude at 
the inner boundary; (M), (u) - mean mass loss rate and outflow veloc- 
ity at the outer boundary; (f c ) - mean degree of condensation of carbon 
into dust at the outer boundary. The radial coordinates in this work are 
plotted in units of the corresponding stellar radii R+ of the hydrostatic 
initial models, calculated from their luminosities and temperatures 
T+ (as given in the table) via the relation L* = 4nR\ crT^. 



2. Model atmospheres and spectral synthesis 

2. 1 . General remarks 

Modelling the cool and very extended atmospheres of evolved 
AGB stars remains challenging due to the intricate interaction of 
different complicated phenomena (convection, pulsation, radi- 
ation, molecular and dust formation/absorption, acceleration of 
winds). Dynamic model atmospheres are constructed to simulate 
and understand the physical processes (e.g. mass loss) occuring 
in the outer layers of AGB stars. In particular they are needed 
if one is interested in reproducing the complex and temporally 
varying atmospheric structures that form the basis for radiative 
transfer calculations, which allows us to simulate observational 
results (spectra, photometry, etc.). 

For our line profile modelling we used dynamic model atmo- 
spheres as described in detail by Hofner et al. (2003; DMA3), 
Gautschy-Loidl et al. (120041 DMA4), N05 or Papers I+II. These 
models represent the scenario of pulsation-enhanced dust-driven 
winds (cf. Sects. 4.7+4.8 of GH04). They provide a consistent 
and realistic description from the deep and dust-free photo- 
sphere (dominated by the pulsation of the stellar interior) out 
to the dust-forming layers and beyond to the stellar wind region 
at the inner circumstellar envelope (characterised by the cool, 
steady outflow). This is accomplished by a combined and self- 
consistent solution of hydrodynamics, frequency-dependent ra- 
diative transfer and a detailed time-dependent treatment of dust 
formation and evolution. 

As a result of deep-reaching convection (dredge-up), nucleo- 
synthesis products can be mixed up from the stellar interior of 
AGB stars, resulting in a metamorphosis of the molecular chem- 
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0.5 1 1.5 



phase bol 

Fig. 1. Movement of mass shells with time at different depths of 
the dust-free model W, which exhibits no mass loss and shows a 
strictly periodic behaviour for all layers. Note the different scales 
on the radius-ordinates compared to Fig. [2] The shown trajecto- 
ries represent the points of the adaptive grid at a selected instance 
of time (higher density of points at the location of shocks) and 
their evolution with time (cf. the caption of Fig. 2 in DM A3 for 
further explanations). 



istry of the whole atmosphere (e.g. Busso et al. 119991 Herwig 
2005). The most important product of all the nuclearly processed 
material mixed up is carbon 12 C. As a consequence, the stars 
can turn from oxygen-rich (C/(D<1) to carbon-rich (C/0>1) dur- 
ing the late AGB phase. The resulting so-called carbon stars 
(C stars) can be found close to the tip of the AGB in observed 
colour-magnitude diagrams (e.g. Nowotny et al. 2003a). The 
drastic change of the atmospheric chemical composition is not 
only relevant for the observable spectral type of the star (chang- 
ing from M to C), it is also crucial for the formation of circum- 
stellar dust. 

Observational studies revealed a rich mineralogy in the dusty 
outflows of O-rich objects (e.g. Molster & Waters 2003 and ref- 
erences therein). Unfortunately, the dust formation process is 
theoretically not fully understood. Many physical and chemical 
details of the process are still not clear, which would be neces- 
sary for a fully consistent numerical treatment. In addition to 
the lack of a grain formation theory, there is an ongoing de- 
bate concerning the underlying physics of the driving mecha- 
nism (Woitke |2006bjtj2007l H6fner|2003 Hofner & Andersen 
120071 ). A potential solution was recently suggested by Hofner 
(120081 . 

Things are quite different for C-rich stars, where we find a 
rather simple composition of the circumstellar dust. A very lim- 
ited variety of dust species were identified by their spectral fea- 
tures (e.g. Molster & Waters 2003), as for example SiC (promi- 
nent feature at »Tl//m) or MgS (broad emission band around 
30yum). The most important species is amorphous carbon dust, 




-0.5 0.5 1 1.5 



phase bol 

Fig. 2. Same plot as Fig. Q] for model M, representing the sce- 
nario of a pulsation-enhanced dust-driven wind. The plot illus- 
trates the different regions within the atmosphere of a typical 
mass-losing LPV. The innermost, dust-free layers below ~2/?* 
are subject to strictly regular motions caused by the pulsating in- 
terior (shock fronts). The dust-forming region (colour-coded is 
the degree of dust condensation / c ) at &2-3 R* where the stellar 
wind is triggered represents dynamically a transition region with 
moderate velocities, not necessarily periodic. A continuous out- 
flow is found from »4 /?+ outwards, where the dust-driven wind 
is decisive from the dynamic point of view. 



though. Not producing any distinctive spectral feature, grains of 
carbon dust represent the dominating condensate and play a cru- 
cial role from the dynamic point of view (mass loss process). 
Amorphous carbon fulfills the relevant criteria for a catalyst of 
dust-driven winds: (i) made up of abundant elements, (ii) simple 
and efficient formation process, (iii) refractory, i.e. stable at high 
temperatures (aT500K), (iv) large radiative cross section around 
lyum in order to absorb momentum. Moreover, the formation 
and evolution of amorphous carbon dust grains can be treated 
numerically in a consistent way by using moment equations as 
described in Gail & Sedlmayr (1 19881 1 and Gauger et al. dl 9901 1. 
As a consequence of this, atmospheric models for C-rich AGB 
stars with dusty outflows were quite successfully calculated and 
applied in the past by different groups (see the overviews given 
by Woitke 2003 or Hofner et al. [2005 ), among them the models 
used in this work. 



2.2. Models for atmospheres and winds 

Following the previous remarks, we concentrated in this study 
on model atmospheres for C-type LPVs for our line profile mod- 
elling, because these models contain a more consistent prescrip- 
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8.0 



Dynamic model atmosphere (C-rich) 



° 7.0 



0.99 



1 — i — i — i — r 



Pulsating AGB star - no dust/mass loss 




0.9 1.0 



1.1 1.2 1.3 1.4 1.5 
R [ R J 



1.6 



S 2 



~i — 1 — i — 1 — i — 1 — i — 1 — i — 1 — i — 1 — r 



hydrostatic initial model 

+ several dynamic phases <p bl 




0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 
R [R H< = 355R ] 




1.1 1.2 1.3 

R [RJ 



1 .6 



Fig. 3. Characteristic properties of model W. Plotted in panel (a) is the bolometric lightcurve resulting from the variable inner 
boundary (piston), diamonds mark instances of time for which snapshots of the atmospheric structure were stored by the radiation- 
hydrodynamics code. Denoted are selected phases for which line profiles are presented in Figf20] Other panels: Atmospheric 
structures of the initial hydrostatic model (thick black line) and selected phases 0boi of the dynamic calculation (colour-coded in the 
same way as the phase labels of panel a). Plotted are gas temperature (b), gas density (c), and gas velocities (d). Note that this is a 
dust-free model and the degree of condensation of carbon into dust is f c = for every depth point at each instance of time. 



tion of dust formation^] However, the velocity effects on line 
profiles are of general relevance, and a comparison with obser- 
vational results of stars with other spectral types (e.g. ^Cyg) 
should be justified. This is especially the case for CO lines 
because of the characteristic properties of this molecule (e.g. 
Sect. 2.2 in Paper I). TableQJists the parameters and the resulting 
wind properties of the dynamic models used in this work. 

The starting point for the calculation is a hydrostatic initial 
model, which is very similar to classical model atmospheres 
(Fig. 1 in DMA3), as for example those calculated with the 
MARCS-Code (Gustafsson et al. 2008). These initial models are 
characterised by a set of parameters as listed in the first part of 
TableQ] They are dust-free and rather compact in comparison 
with fully developped dynamical structures (e.g. Fig.|4|i. The ef- 
fects of pulsation of the stellar interior are then simulated by a 
variable inner boundary 



Au„P (2tt 
Rm(f)=«fa(0) + =^sinlyf 



(1) 



1 First M-type models were presented by Jeong et al. (2003 ). 
However, the development of a wind driven by radiation pressure on 
dust grains of O-rich species in these models possibly results from an 
overestimation of dust opacities in the course of grey radiative transfer 
as applied by these specific models (Woitke 2006b] H6fner[2007 ). 



This so-called piston moves sinusoidally with a period P and a 
velocity amplitude Au p . Assuming a constant radiative flux at 
the inner boundary, the luminosity Li„(t) varies like oc /??(?) for 
the models presented in DMA3. The assumptions for the inner 
boundary were slightly adapted for the latest generation of mod- 
els as discussed in DMA4: an additional free parameter /l was 
introduced so that the luminosity at the inner boundary varies 
like 



Lm(t) - Lp 

Lo 



R (?) - K 

m v ' 



2-| 



R 2 



(2) 



where Lo and Ro are the values at the inner boundary of the hy- 
drostatic initial model (Lo=L+, by definition, and /l=1 corre- 
sponds to the previously used case of constant flux at the inner 
boundary)Q The luminosity variation amplitude of the model 
can be adjusted thereby independent of the mechanical energy 



2 Note that the average L over a period is slightly higher in a dy- 
namical model than in the corresponding hydrostatic initial model due 
to the quadratic dependence of L m {t) on R m {t). The effect is small for 
typical pulsation amplitudes,though («0. 1% for model W, <2% for 
model M). Note also that R =R m {Q) is not equal to R+ (calculated as 
described in Tab.fTJ, but somewhat smaller, as the inner boundary is lo- 
cated below the photopshere (e.g. for model S: i? =2.832157-10 13 cm, 
R* =3.42803- 10 13 cm). 
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Dynamic model atmosphere (C-rich) 




Pulsating AGB star (Mira) with mass loss 

-srrT 



2 3 4 

R [R=412R ] 




cn 




Fig. 4. Characteristic properties of model M. Radial structures of the initial hydrostatic model (thick black) and selected phases 
of the dynamic model atmosphere during one pulsation cycle (4>b o \=0.0-1.0). Plotted are gas temperatures (a), gas densities (b), 
gas velocities (c), and degrees of condensation of carbon into dust (d). Panel b demonstrates the larger extension compared to the 
hydrostatic initial model and the local density variations. 



input by the piston. This allows us to tune the bolometric am- 
plitudes Amboi to resemble more closely the values derived from 
observational studies. Figure^ shows the variability of the lu- 
minosity L of one dynamic model as an example. 

As described in detail in DMA3, the models can be divided 
into two sub-groups according to their dynamical behaviour: 

> pulsating model atmospheres where no dust forms 

> models developing pulsation-enhanced dust-driven winds. 

Model W (directly taken from DMA3 where it is denoted 
by 170t28cl4u2) represents a typical example for the first type 
of dynamical models. FigureQ] shows the temporal evolution of 
mass layers in different atmospheric depths for such a dust-free, 
pulsating atmosphere. A mild pulsation is simulated by the small 
amplitude of the piston in this case. However, the matter is not 
ejected far enough so that dust formation is made possible by 
the low temperatures. There is no dust-driven outflow, and pul- 
sation is the sole reason for atmospheric dynamics. The outer 
boundary follows the movement of the upper atmospheric lay- 
ers in the case of models without mass loss. A shock wave arises 
once in every pulsation period, causing the layers to follow more 
or less ballistic trajectories. The resulting completely periodic 
behaviour can be recognised in Fig.Q] The atmospheric struc- 
ture for different phases throughout a pulsational period is shown 
in Fig. [3] Varying around the hydrostatic configuration, the dy- 
namic model does not resemble it at any point. 



The atmospheres of LPVs become qualitatively different 
with the occurence of dust and the development of a stellar wind. 
The other models listed in TableQ] namely models S and M, 
represent examples for this second type of dynamic model at- 
mospheres. Again, shock waves are triggered by the pulsating 
stellar interior and propagate outwards. The difference arises be- 
cause efficient dust formation can take place in the wake of the 
shock waves (post-shock regions with strongly enhanced den- 
sities at low temperatures). Radiation pressure on the newly 
formed dust grains results in an outflow of the outer atmo- 
spheric layers. This behaviour is demonstrated in Fig. [2] which 
shows the characteristic pattern of the moving mass layers with 
model M as an example. While the models without mass loss 
stay rather compact compared to the hydrostatic initial model 
(Fig.O, the models that develop a wind are inflated by it and be- 
come much more extended than the corresponding initial model. 
Numerically, a transmitting outer boundary - allowing outflow; 
fixed at 20-307?+ - is used for the latter type of dynamic model 
atmospheres. Spatial structures of model M are shown in Fig. [4] 
demonstrating the strong influence of dust formation on the at- 
mospheric extension. The models S and M exhibit quite a mod- 
erate dust formation process, which leads to a smooth transistion 
of the velocity field from the pulsating inner layers to the steady 
outflow of the outer layers (see the velocity structure plots in 
Fig. 9 of Paper I and Fig. [4] in this work). In contrast, the model 
applied in Sect. 6.1 of Paper II is a representative of a group of 
models with a more extreme dust formation. Not every emerg- 
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ing shock wave leads to the formation of dust grains (cf. Fig. 2 
in DMA3) and the velocity field in the dust-forming region may 
look very different for similar phases of different pulsation peri- 
ods. Pronounced dust shells arise from time to time and propa- 
gate outwards (see the structure plot in Fig. 10 of Paper II). 

We refer to Hofner et al. (2003 ) for more details about the 
numerical methods and the atmospheric models. 

2.3. The specific models used 

Model S was already used extensively in the previous Papers I+II 
to study line profile variations. Its parameters (listed in Table. [TJ 
were chosen to resemble the Mira S Cep, as this is the only 
C-type Mira with an extensive time series of high-resolution 
spectroscopy and derived RVs (Fig. [XT}. Table 1 of Paper II lists 
properties for this object as compiled from the literature for com- 
parison. Model S should not be taken as a specific fit for S Cep, 
though. A discussion of the difficulties when relating dynamic 
model atmospheres to certain objects can be found in Sect. 3 
of Paper II (or Sect. 2.1.4 in N05). However, there is evidence 
that the model reproduces the outer layers of this star reasonably 
well. There are some properties listed in the mentioned tables 
which can be compared directly and agree to some extent (e.g. 
L, P, M, M eX p)- In addition, low-resolution synthetic spectra in 
the visual and IR computed on the basis of model S resemble ob- 
served spectra (ISO, KAO) of S Cep fairly well as it was shown 
by Gautschy-Loidl et al. (2004; their Sect. 5.1). 

Model M was the (so far) last model of a small parame- 
ter study we carried out subsequent to Papers I+II. The inten- 
tion was not to find a model fitting a certain target of observa- 
tions better, but to change the model parameters in order to re- 
produce one particular observational aspect: namely the veloc- 
ity variations in the inner, dust-free photosphere resulting in a 
rather uniform RV curve for CO Av = 3 lines in spectra of Miras 
as discussed in Sect. l6.1l and shown in Fig. [14] The comparison 
of the global velocity field of this model M with observational 
results in Sect. 16. 21 is still done on a qualitative basis, though. 
Confronting RV measurements (Fig.fTTI) with the corresponding 
synthetic values (Fig.[T8l) provides some general information of 
how realistic the velocity structures of the model are. However, a 
direct and quantitative comparison of stellar and model parame- 
ters is not feasible at the moment due to limitations on the obser- 
vational side (very small number of stars observed extensively, 
often only rough estimates for properties of targets) as well as 
on the modelling side (only C-rich models, laborious process to 
get a RV diagramm as Fig.[T8l) as discussed in Sect. 16. 21 

While investigating synthetic CO Av = 3 line profile varia- 
tions based on a few dynamical models available at that time, 
we found model W reproducing the observed behaviour of 
SRVs CSect. 17. It with similarities to WHya (increased ARV). 
In Sect. 17.31 we will make a comparison of observational and 
modelling results only for this very selected aspect. Relating 
model W (C-rich chemistry) and the M-type LPV W Hya to con- 
strain the parameters of this star is even less possible than for the 
before mentioned case of model M, and was not intended. 



2.4. Calculating synthetic line profiles 

For the spectral synthesis we followed the numerical approach 
as described in detail in N05 (Sect. 2.2) and also in Papers I+II. 
The modelling procedure of dynamic model atmospheres (see 
DMA3) yields some immediate results, like mass loss rates M, 
terminal velocities of the winds Voo, or degrees of dust condensa- 



tion in the outflows f. In addition, it provides snapshots of the 
time-dependent atmospheric structure (p, T, p, u, etc.) at several 
instances of time (e.g. Figs.[3}t©. These represent the starting 
point for the aspired spectral synthesis, accomplished in a two- 
step process as described below. 

The first step, calculating opacities based on a given atmo- 
spheric structure, was accomplished with the COMA code, a de- 
scription of which can be found in Aringer (2000), Gautschy- 
Loidl (2001), or N05. Informations on recent updates and the 
latest version can be found in Gorfer (2005), Lederer & Aringer 
(2009), Aringer et al. (2009). Element abundances for the spec- 
tral synthesis were used in consistency with the hydrodynamic 
models of solar composition. We adopted the values from 
Anders & Grevesse ( 119891 ), except for C, N and O where we 
took the data from Grevesse & Sauval ( 1994). This agrees with 
our previous work (e.g. Aringer et al. 1999, Aringer et al. 2009) 
and results in Z « 0.02. Subsequently, the carbon abundance 
was increased according to the C/O of the models. Abundances 
and ionisations of various atoms and formed molecules for all 
layers of the atmospheric model were calculated with equilib- 
rium chemistry routines (for a detailed discussion and an ex- 
tensive list of references we refer to Lederer & Aringer 2009 ). 
The depletion of carbon in the gas phase due to consumption by 
dust grain formation is also taken into account by the COMA 
code. Examples for the resulting partial pressures can be found 
in Fig. [5] Subsequently, opacities for every radial depth point 
and the chosen wavelength grid were computed. Several opac- 
ity sources were considered, the most important one being the 
molecules for which the line profiles are to be studied. Their 
contribution to the opacities were computed by using line lists 
and under the following assumptions: (i) conditions of LTE, (ii) 
a microturbulence velocity of ^=2.5kms~', and (iii) line shapes 
described by Doppler profiles. Assuming LTE conditions, level 
populations can be computed from Boltzmann distributions at 
the corresponding gas temperature T . The opacity k y at a given 
frequency v for a certain transition from state m to state n can 
then be written as 



Nne 1 gf 



-Eo/kT 



{l-e- hVolkT )<p(v), 



»>,<■ Q(T)' ' {J) 

with the number density of relevant particles N, the charge e and 
mass m e of the electron, the speed of light c, the Boltzmann con- 
stant k, and the energy hvo of the respective radiation. The par- 
tition function Q(T) is the weighted sum of all possible states. 
Eq represents the excitation energy of the level m (from ground 
state) and gf is the product of the statistical weight g( m ) of the 
level times the oscillator strength /(,„,„) of the transition. Line 
lists usually contain frequencies vo (or in practice wavenum- 
bers), excitation energies Eo, and gf values together with in- 
formations for line identification. In order to reproduce the line 
shapes in a realistic way, a broadening function ip(v) for the line 
profile was introduced with 



f 

Jo 



tp(v) dv — 1 



(4) 



Only the effects of thermal broadening (first term in Eq.[6]l and 
the non-thermal contribution of microturbulent velocities (sec- 
ond term in Eq.|6]l are taken into account by COMA, whereas 
other effects (e.g. natural and pressure broadening, macroturbu- 
lence) are neglected. We refer to Sect. 2.2.2 of N05 for details. 
The resulting Doppler profiles can be described by a (Gaussian) 
broadening function 



<p(y) 



1 



(5) 
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with a Doppler width A v given by 



A V 

A v = — 



V 



2HT 



+ e 



(6) 



where % is the gas constant, fi the molecular weight, and £ the 
microturbulent velocity. An updated set of references for the 
line lists of all molecular species used by the current version 
of COMA can be found in Lederer & Aringer (20®. In this 
work we made use of the list of Goorvitch & Chackerian ( 1994 ) 
for CO and the CN line list of J0rgensen & Larsson ( 1990}. We 
utilised the same molecular lines as in Papers I+II, their respec- 
tive properties are summarised in Table[2] In addition, molecules 
affecting the spectra pseudo-continuously (mainly C2H2) were 
considered by an opacity of constant value for the considered 
spectral range (cf. Sect. 5.1. in Paperl or N05). Furthermore, 
continuum absorption coefficients are determined by COMA 
(Lederer & Aringer [20091 1 as well as the opacity due to dust 
grains of amorphous carbon. For the latter we used the data of 
Rouleau & Martin ( 199 1] set AC) and computed the resulting 
dust absorptior0 under the assumption of the small particle limit 
of the Mie theory (grain sizes much smaller than relevant wave- 
lengths -v» opacities proportional to the total amount of con- 
densed material, but independent of grain size distribution; cf. 
DMA3, Hofner et al. [T9921 

In the second step of the spectral synthesis, the previously 
calculated data array of opacities k v (A,t) for all depth and wave- 
lengths points was utilised to solve the radiative transfer (RT). 
As the thickness of the line-forming region in AGB atmospheres 
is large compared to the stellar radii, it is necessary to treat 
the RT in spherical geometry. In addition, the complex veloc- 
ity fields (pattern of outflow and infall as for example shown in 
Fig® in AGB atmospheres severely affect the line shapes in the 
resulting spectra (observed or synthetic, e.g. Fig.[T2l. and it is 
essential to include the influence of relative macroscopic veloci- 
ties in this step of the spectral synthesis. Thus, a code for solving 
spherical RT, which takes into account velocity effects, is used 
to model line profiles and their variations. The RT code used 
in this work (Windsteig 1998) follows the numerical algorithm 
described in Yorke (fT9§8i l. 

It is necessary to choose spectral resolutions that are high 
enough to sample individual spectral lines with a sufficient num- 
ber of wavelength points. This is especially important for syn- 
thesising the often quite complex line profiles for stars with 
pronounced atmospheric dynamics, like those that are the topic 
of this work. Therefore, all spectra are calculated with an ex- 
tremely high resolution of /?=/l/A/l=300000 and were then re- 
binned to /?=70000 for comparison with observed FTS spec- 
tra. Furthermore, the synthetic spectra shown below were nor- 
malised relative to a computation with only the continuous opac- 
ity taken into account (F/F cont ). 

Our aim is to infer information about atmospheric veloc- 
ity fields from shifts in the wavelength of (components of) 
spectral lines. For this purpose, RVs were calculated by using 
the rest wavelength of the respective line and the formula for 
Doppler shift. For an easy comparison of our modelling results 
with observations, we adopted the naming convention for ve- 
locities of observational studies (positive for material moving 
away from the observer and negative for matter moving towards 
the observer). Thus, outflow from the star results in blue-shifted 
lines and negative RVs, while infalling matter revealed by red- 
shifted lines leads to positive RVs. In general, the RVs resulting 

3 Scattering was not taken into account, neither for the dynamic mod- 
els nor for the spectral synthesis. 



Table 2. Different molecular features and the properties of the 
specific lines used for the line profile modelling here (and in 
Papers I+II). 





Specific 


line chosen for modelling 


Line type 


designation cr/wn [cm '] 


A [/mi] 


CN Av = -2 red 


1-3 Q 2 4.5 


4871.3400 


2.0528 


CO Av = 3 


5-2 P30 


6033.8967 


1.6573 


CO Av = 2 high-exc. 


2-0 R82 


4321.2240 


2.3142 


CO Av = 2 low-exc. 


2-0 R19 


4322.0657 


2.3137 


CO Av = 1 


1-0 Rl 


2150.8560 


4.6493 



Notes: The lines were chosen according to the criteria given in Sect. 5.1 
of Paper I. The values given here correspond to the rest wavelengths 
or central frequencies v as used in Eq. Throughout this work we 
will use only a low-excitation first overtone CO line (i.e. CO 2-0 R19) 
denoted by CO Av = 2, while the high-excitation line CO 2-0 R82 is 
merely listed for completeness as it was used in Papers I+II. This dis- 
tinction is necessary as the two types of lines (at least the ones used 
by HHR82 and similar studies) behave differently in kinematics studies 
(see Fig.ll7l>. 



from observations were combined into one composite lightcycle 
(</>v=0. 0-1.0) and then plotted repeatedly for better illustration 
(e.g. Fig.fl~4b. For the modelling we computed spectra and de- 
rived RVs throughout one pulsation period ((f)\, o \-0.0-1.0) and 
replicate the values beyond this interval (e.g. Fig.fTBIl. 

Following the convention of Papers I+II, bolometric phases 
0boi within the lightcycle in luminosity (Fig. [3^) will be used 
throughout this work to characterise the modelling results (at- 
mospheric structures, synthetic spectra, derived RVs, etc.) with 
numbers written in italics for a clear distinction from the visual 
phases V , which are usually used to denote observational results 
(with (p v -Q corresponding to phases of maximum light in the vi- 
sual). For a discussion of the relation between the two types of 
phase informations, namely 0boi and </> v , we refer to AppendixlAl 

Throughout this work (e.g. the lower panel of Fig.|6]and all 
similar plots in the following), the radial optical depth r v (r) is 
computed by radially integrating inward^ for a given wave- 
length point 



Ty(r) 



cont(r') + fv,lines('-') + K v,dust(r')] p(r') d/ 



(7) 



and is then plotted with dotted lines. For the opacities K v (r), 
all relevant sources (continuous, molecular/atomic, dust) are in- 
cluded. If velocity effects are taken into account, optical depths 
are calculated with the opacities at all depth points Doppler- 
shifted according to the corresponding gas velocity there and 
then plotted with solid lines. The optical depth of 7" y «l provides 
clues on the approximate location of atmospheric layers where 
radiation of a certain frequency v (or the corresponding wave- 
length A) originates. 



3. Line formation in the extended atmospheres of 
evolved red giant stars 

This section is concerned with the fact that different molecu- 
lar spectral features visible in spectra of AGB stars originate at 
different atmospheric depths and the examination of this with 
numerical methods. 



4 corresponding to the central beam with impact parameter p=0 in 
the framework of spherical RT (e.g. Fig. 1 of Nordlund ll984t 
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Typical main sequence stars exhibit relatively compact pho- 
tospheres, and the radiation in different wavelengths origi- 
nates at roughly the same geometrical atmospheric depths with 
a rather well defined temperature (e.g. Sect. 4.3.1 in GH04). 
For example, the relative thickness of the flux-forming re- 
gion where the spectrum is produced in our Sun amounts to 
A^phot,o/^o ~ 4-1CT 4 (cf. N05). The often used parameter effec- 
tive temperature T e ff represents a mean temperature in the layers 
of the thin photosphere, from where almost all photons can es- 
cape. 

In contrast, estimating the atmospheric extensions AR atm + in 
the case of AGB stars is much more difficult. The outer bound- 
aries of these objects are hard to define (as discussed by GH04 in 
their Sect. 4.1) due to the shallow density gradients related to the 
low surface gravities g+, intensified by the dynamic processes 
of pulsation and mass loss. GH04 give estimates of 0.1-0.5 for 
the parameter A/? atmjlk /l? + . For the dynamic models used here 
('Sect. l2.2t . one could for example consider the region of dust 
formation at ~2/?+, where the onset of the stellar wind takes 
place (cf. FigUJi as an outer boundary of the atmosphere. This 
results in AR atm .+//?* m 1 . Whatever value one adopts, it is clear 
that the thickness of an AGB atmosphere is of the same order as 
the stellar radius, which is quite different from normal main se- 
quence stars. With these extremely extended atmospheres, the 
wavelength-dependent geometrical radius can vary by a few 
100 R Q (equivalent to a few AU) or more. Also the related gas 
temperatures can cover a wide range from «5000K down to a 
few hundred Kelvin (N05). 

In general, molecular lines do not originate in a well-defined 
and narrow region, but are formed over a more or less wide range 
in radius. Nevertheless, there have been attempts to constrain 
the approximative line forming region in radius and temperature 
for different molecular features with similar dynamical model 
atmospheres as used here. The results - which should be re- 
garded as rough estimates rather than definite borders, though 
- are compiled e.g. in Table 1.2 of N05. There are only a few 
molecular bands (e.g. CN or C2), which are solely formed in the 
deep, warm layers of the photosphere. Most of the strong fea- 
tures (e.g. H2O, HCN, C2H2) originate in the cool, upper layers 
of the atmosphere because of the small dissociation energies of 
the (polyatomic) molecules in combination with the high gf val- 
ues. It is interesting that there are examples for bands formed in 
a relatively narrow region, as e.g. the ones of C3, which is due 
to the interplay between molecular formation and dissociation 
probabilities (Gautschy-Loidl 2001 Fig.[5]l. The special role of 
CO and the corresponding features will be discussed in detail 
below. To sum this up, the emerging spectrum (across the whole 
spectral range from the visual to the IR) of an AGB star is formed 
over a large radial range with diverse physical conditions, which 
is quite different from the scenario for the Sun as sketched above. 

Moreover, not only spectral features of different molecules 
originate at varying radial atmospheric depths, but also different 
lines of the same molecular species can show this effect. This 
is especially pronounced in the case of CO lines and will be in- 
vestigated below on the basis of one (arbitrarily chosen) phase 
of model S. The upper panel of Fig. [5] shows the corresponding 
atmospheric structureQat the phase (p\, o i=0.0. 



3 Note: The dynamical models reach into geometrical depths with 
temperatures of about 4-5000 K (Fig.Q. Modelling even deeper zones 
would (besides numerical difficulties) mean to include convection and 
the driving mechanism of pulsation, which has not been done so far. 
The resulting overall continuous optical depth is too low (<1) at cer- 
tain wavelengths, but this is not the case for the wavelengths of strong 




-10 I 

-12 
-14 o 
-16c 



.7 1 2 3 5 7 10 

radius R [R >[ =493R S ] 

Fig. 5. Upper panel: Atmospheric structure of model S for the 
phase (pb o \-0.0. Gas velocities overplotted for illustration 

purposes (shock fronts). Compare Figs. 1 and 9 in Paper I to get 
an idea of the temporal variations and also for the scale of the gas 
velocities. The dotted line marks the border between the actual 
model and the extension towards the interior for the sake of opti- 
cal depth (see text). Lower panel: The total gas pressure together 
with the corresponding partial pressures of selected molecules 
for this phase if chemical equilibrium is assumed. Shown is only 
a subset of all molecular species considered by COMA, which 
are CO, CN, C 2 , CH, C 3 , C 2 H 2 , and HCN. 

There are a few factors that are decisive for the line intensi- 
ties and for the locations (r, r eas ) of the line forming regions of 
molecular features in AGB atmospheres: 



111. 



the atmospheric structure itself (r-T-p) 
the relative abundance (i.e. partial pressure) of the consid- 
ered molecular species at certain depth points within the at- 
mosphere (formation/depletion) 

the properties of the respective line (excitation energy E exc , 
strength of the transition gf) 

possibly other opacity sources (local continuous opac- 
ity, other molecular lines, dust opacity, pseudo-continuous 
molecular contributions) 

influence due to the relative macroscopic velocity fields. 



Several studies demonstrated in the past that the atmo- 
spheric structure (— >i.) has substantial impact on spectral fea- 
tures. Compared with hydrostatic model atmospheres, dynamic 
models have enhanced densities in cool upper layers (Fig.|4j». 



absorption features, as illustrated in Nowotny et al. (2003b). Thus, a 
linear extrapolation in T-p to deeper layers (r gas «10 4 K, t«100) had 
to be applied for the spectral synthesis presented here. Velocities were 
considered to be zero in the extension. This assumption is not problem- 
atic since these layers are only relevant for continuum optical depths 
and have no influence on the line profiles as molecular lines are formed 
further out. 
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This results in a change of molecular abundance and decreased 
or increased intensities of molecular features (as outlined in 
Sects. 4.7.5 and 4.8.3.2 of GH04). Examples for this effect would 
be the SiO bands at 4yum (Fig. 5 of Aringer [T999l l. the H 2 bands 
around 3/im (Fig. 9 in DMA3), or the combined feature of C2H2 
and CN at * 14/xm (Sect. 4 of DMA4). 

The lower panel of Fig. [5] demonstrates the differences in 
relative abundance (— > ii.) for various molecules throughout the 
whole atmosphere. Plotted are the partial pressures for differ- 
ent molecular species resulting from the evaluation under the 
assumption of chemical equilibrium by the COMA code. CN 
is a typical representative for species which form in deep pho- 
tospheric layers below »T.5/?*. Values of 2400-5000K for the 
gas temperatures in the line forming regions were assigned to 
this species by Gautschy-Loidl (120011 1. In contrast, C2H2 shows 
a quite different behaviour. As they are very sensitive to temper- 
ature (dissociation), these molecules can be found from w 1 .2 
on outwards, Gautschy-Loidl lists temperatures of < 2400 K. C3 
represents a fairly extreme example. Its formation requires at the 
same time high densities for high enough collision probabili- 
ties with low temperatures in order not to dissociate. This is met 
only in a rather narrow region around temperatures of 2000 K. 
The density variations due to the propagating shock wave (for 
this phase at a radius slightly larger than 1R+) also are of im- 
portance, as can be seen in the partial pressures. From ^0.8/?* 
inwards the abundances of all species decrease as molecules are 
dissociated for temperatures higher than 5000 K. 

Figure|5] illustrates also that CO is an outstanding molecular 
species concerning its relative abundance. Due to its high dis- 
sociation energy (11.1 eV), CO is stable and the most abundant 
species (of those shown; not included in this figure is e.g. the by 
far most abundant molecular species H2) at all depth points of 
the atmospheric model. For some molecules - like CN or even 
more C3 - the limited range with significant large partial pres- 
sures sets constraints on the line forming region. However, CO 
is present and abundant across the whole atmosphere, and the 
properties of individual spectral lines (— >iii.), namely the ex- 
citation energy E exc and the strength of the transition gf, be- 
come important. The excitation energy Eq of the lower level of 
a certain transition m — > n determines at which temperatures (i.e. 
atmospheric depths) the levels are populated so that an absorp- 
tion can occur at all. In addition, the strength of the transition 
(g/-value) determines the intensity of the line and thus also the 
atmospheric depth were r y «l is reached. Compare Sect. 12. 41 and 
Eq. (01 in this context. Also the relative abundance is present in 
this equation by the quantity N. 

To study these effects for the lines chosen for our line profile 
modelling (cf. Table[2] Fig. [8]), we used the atmospheric struc- 
ture shown in Fig. [5] and calculated different quantities. Note 
that velocity effects were neglected for all subsequent compu- 
tations of this subsection. The results are plotted in Fig. [6] The 
upper panel of this plot shows the molecular absorption (ab- 
sorption per mol of the material of the corresponding molecu- 
lar species), illustrating how one individual molecule of a given 
species would absorb radiation of the wavelength of the corre- 
sponding line at all atmospheric depth points. Level populations 
are computed by Boltzmann distributions for the respective tem- 
perature (LTE). The increasing gf values from second overtone 
(Av = 3) to first overtone (Av = 2) and to fundamental (Av= 1) 
CO lines (cf. Fig. [8]) are reflected in the increasing molecular 
absorption. The middle panel of Fig. [6] shows the actual absorp- 
tion coefficient (absorption per gram of the whole stellar ma- 
terial) due to the respective molecular line, where the relative 
abundance of the species comes in. This can be recognised by 




Fig. 6. Radial distribution of absorption per mol of the molecu- 
lar material (upper panel), absorption coefficient (centre panel) 
and the optical depths (see text for detailed discussion) for the 
molecular lines chosen for the line profile modelling as listed 
in Table|2] calculated on the basis of the atmospheric structure 
shown in the upper panel of Fig. [5] For the absorptions (upper 
+ middle panel) only the contribution of the respective lines are 
plotted, while all relevant opacity sources (cf. Eq.[7]) are taken 
into account for the optical depths (lower panel). Note that the 
effects of gas velocities (Doppler-shifts) are not taken into ac- 
count for the computation of the optical depths. 



the decrease for all k v values with higher temperatures or the 
steep decrease for the CN line below 2000 K where the partial 
pressure of CN severely drops (Fig.|5]l. Plotted in the middle 
panel is the opacity due to the CO or CN lines only. In com- 
bination with all other contributions (continuous, other molec- 
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Fig. 7. Same as bottom panel of Fig. [6] the individual contribu- 
tions of different opacity sources (continuous, molecular line, 
dust) to the radial optical depth (cf. Eq.[7|l for two of the CO 
lines used in this study (Tab.[2]i are shown in addition. 



ular species, dust; — >iv.), these opacities serve as input for the 
RT or can be used to calculate radial optical depths as shown in 
the lower panel of Fig. [6] The quantities k v and r v for the cen- 
tral/rest wavelength of the respective lines (cf. line profiles in 
Fig.[T2lplotted grey) are shown here, no velocity effects are taken 
into account. It can clearly be seen that t v «T is reached in a wide 
temperature range between «300 K and «3000 K, corresponding 
to A/telO-20fl**23-46At/ in this case. While for the CN line 
the limited region of high enough partial pressure and the line 
parameters are crucial for the line formation, only the latter is 
relevant for CO lines. Different rotation-vibration band systems 
of CO originate in quite separated regions within the stellar at- 
mosphere. Therefore and because of other important effects (cf. 
Sect. 2.2 of Paper I), the spectral lines of CO play a major role 
in the research on AGB atmospheres, especially for kinematic 
studies. 

The contributions of different opacity sources to the result- 
ing radial optical depths are illustrated in Fig. [7] for the model 
atmosphere used (Fig. [5]). The upper panel of this plot shows 
that the strong fundamental CO line is the dominant source at 
its central/rest wavelength over almost the whole extended at- 
mosphere. Only in the hot layers below the photosphere the con- 
tinuous opacity takes over, while the dust plays a neglectable 
role in general. The panel below shows the same quantities for 
the chosen second overtone CO line. Because of its properties 
(lower gf value and higher £ e xc — > lower line intensity com- 
pared to the Av = 1 line; cf. the spectrum in Fig. [8]), this molecu- 
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Fig. 8. Characteristic properties excitation energy and gf value 
of several lines of the CO Av = 3 / 2/ 1 bands are shown in the up- 
per and middle panel. Plotted in each case are lines of the same 
rotational transition but from different vibrational bands (mean- 
ing for example 1-0 Rl, 2-1 Rl, 3-2 Rl, ... , 21-20 Rl). The 
strong variations in E exc and gf from one vibrational band to 
another influence the results shown in Fig.|6]and lead to distinct 
different line intensities. This is illustrated in the lower panel 
with a spectrum (frequency times specific luminosity vs. wave- 
length) based on the atmospheric model of Fig.|5]where only the 
respective lines are taken into account. Marked with the same 
colour-code as Fig. [6] are the lines actually used for the line pro- 
file modelling (as listed in Table|2]i. 



lar line provides the main contribution to the total optical depth 
only in the layers of w 2-4000 K. From « 5000 K inwards the 
continuous opacity represents the major source again, although 
to a somewhat lower extent (minimum of absorption due to the 
H~ ion around 1.6pm; e.g. GH04) compared to the former case 
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Fig. 9. Same as bottom panel of Fig.[6]but for the large number 
of CO Av = 3/2/1 lines as shown in Fig. [8] Plotted with thick 
linestyle are the CO lines that were actually used for the line 
profile modelling (Fig[6]). 



(^4.65 pm). On the other hand, the absorption due to grains of 
amorphous carbon dust is stronger in the wavelength region of 
CO Av = 3 lines (cf. Fig. 1 of Andersen et al. 119991) as reflected 
in the dust optical depths. Still, the absorption features are vis- 
ible in the synthetic spectra because of the moderate mass loss 
of the used model S (Tab.[T]i. This changes for higher mass loss 
rates, the intensities of molecular features decrease as the corre- 
sponding line forming regions are hidden by optically thick dust 
shells. 

Figure|8] illustrates that the same type of CO lines (funda- 
mental, first or second overtone), but from different vibrational 
bands can have very different values of E exc and gf. Not only 
leading to variations in line intensities (lower panel), it also has 
a strong impact on the absorption coefficients and optical depths 
as discussed above. This is demonstrated in Fig.|9l where the 
radial optical depths are plotted for the variety of lines from dif- 
ferent vibrational bands used for Fig. [8] Going to higher vibra- 
tional quantum numbers leads to weaker lines and line formation 
regions located at higher temperatures (T gas ), i.e. further inside 
the atmosphere. One may suppose that these differences should 
also be reflected in the dynamic effects (line profile variations, 
derived RVs) if other lines than those used for our line profile 
modelling (Tabled chosen to be similar to observational stud- 
ies) are considered. 

The remaining factor (— >v.) concerning the line formation 
process, namely the important point of the influence of velocity 
fields, will be the main issue in the next sections. Velocity effects 
are most relevant for line profiles in spectra of evolved red giants. 

4. The influence of velocity fields 

The previous section covered the line formation process within 
the extended AGB atmospheres with emphasis on the aspects 
of atmospheric structures, molecular abundances, and properties 
of molecular lines. Now, we will introduce another important 
item relevant for the line formation in the atmospheres of long- 
period variables. Representing the major topic of this work, at- 
mospheric velocity fields with relative macroscopic motions of 
the order of lOkms -1 are decisive for the resulting profiles of 
individual spectral lines. This is illustrated below by numeri- 
cal computations for selected examples of increasing complex- 
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Fig. 10. Synthetic CO Av = 2 line profiles based on the hydro- 
static initial atmospheric structure of the dynamic model W 
(Fig. [3]). Artificial gas velocities constant for each atmospheric 
depth point were assumed, these are listed in the legend on the 
right in units of [kms _I ] and are colour-coded in the same way 
as the spectra. The velocity field has a strong influence on the re- 
sulting line profiles (shapes, intensities, shifts, emission compo- 
nents). Outflow results in blue-shifted absorption profiles, while 
infall leads to red-shifted ones. Values for the conversion factor 
p, as derived from the Doppler-shifted line profiles shown here, 
are listed in Table[3] 

ity. Compare also Mihalas (1978; Sect. 14) in this context for a 
discussion of the theoretical background. 

4.1. The simple case 

The first example is shown in Fig.[l0l representing the most 
simple case conceivable for illustration purposes. For this, we 
used the hydrostatic initial atmospheric structure of model W 
(cf. Table[T]i. Based on this atmospheric structure, synthetic line 
profiles for a selected CO Av = 2 line (Tabled were calculated 
under the assumption of uniform velocity fields, i.e. spherically 
symmetric outflow or infall with one constant radial velocity for 
every atmospheric depth point. This leads to profiles shifted to 
the blue or red, according to the assumed outflow or infall, re- 
spectively. However, due to projection effects the gas velocities 
do not convert directly into observed Doppler-shifts and the mea- 
sured RVs are always smaller. The shapes of the absorption fea- 
tures change obviously, they become broader and more cone-like 
but less deep. This effect is increasing with higher velocities. In 
addition, a weak emission is appearing, Doppler-shifted in the 
opposite way than the absorption component. This model rep- 
resents an extended red giant star with outer layers that are not 
optically thick. Therefore also the surrounding material can sig- 
nificantly contribute to the spectrum, leading to this typical pro- 
file of P Cygni-type when velocity effects are accounted for. 

A correction factor p is needed to connect RVs of Doppler- 
shifted spectral features with actual gas velocities in the corre- 
sponding line forming region by M gas = p • RV. Several studies 
(e.g. Willson et al. 1982 Scholz & Wood 120001 1 investigated in 
the past the size of p, as outlined in Sect. 4.6 of Paper II. We 
also derived values for p for the line profiles shown in Fig.fTUl 
For this purpose, the assumed gas velocities M gas were divided 
by radial velocities RV as measured from the deepest point of 
the line profile. The results, summarised in Table[3] are consis- 
tent with the findings from the literature listed above (and more 
recent results concerning Cepheids by Nardetto et al. 120071 or 
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Table 3. Conversion factors p as derived from the line profiles 
shown in Fig.fTOl 
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Groenewegen 2007). By applying the same method to a hydro- 
static MARCS model atmosphere from Aringer (120001 1 resem- 
bling the star Arcturus (K1.5III), we found values reaching up 
to p ~ 1.5 (see Sect. 3.4 of N05). 

In some cases, the shapes of observed spectral lines may ap- 
pear relatively similar to those shown in Fig.[l0] One example 
is given in the right panel of Fig. [20] with the profiles of lower 
resolution based on model W. In general, however, the velocity 
fields in AGB atmospheres are more complex than uniform ve- 
locity fields and the resulting line profiles show - especially for 
high spectral resolutions - effects beyond shifts in wavelength. 
This will be examined below. The estimation of the conversion 
factor p can also be hindered as illustrated in Sect. [8] 
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4.2. Consistent dynamical AGB star spectra 

The second example is shown in Fig.[TT] where we present a 
synthetic CO Av = 3 line profile based on one selected phase of 
model W ((p\, o \=0.48; part of the results in Fig. 120)1. The velocity 
structure of the chosen phase in the upper panel shows infall for 
all atmospheric points. However, the velocities are quite diverse 
and not even monotonic, especially in the line forming region 
of the second overtone CO lines below ~l/?+ (cf. Sect. l7.4t . 
The smearing of opacity due to varying Doppler-shifts at dif- 
ferent depths leads to a clearly asymmetric shape of the ab- 
sorption line, as shown in the insert of the lower panel. This 
is in addition illustrated by plots of the radial optical depth at 
wavelengths corresponding to the line center (A), the line wing 
(B), and a (pseudo-)continuum point next to the line (C) in the 
lower panel. The bulk of the opacity in the relevant layers is red- 
shifted and the deepest point of the spectral line (B) is located at 
RV- 3.55 kms -1 if velocities are taken into account. The layers 
directly below the emerging shock wave (<0.95 R+) contribute 
most to the wavelengths marked with (A) and (B). In addition, 
there is a strongly enhanced absorption at wavelengths corre- 
sponding to RV~8 km s _1 (C). This is contributed by layers with 
logpg as <-l 1.5 g cirT 3 (cf. Sect. 3.3.1 in N05) slightly outwards 
the step in density at =s0.95,R* and leads to the destorted line 
profile. 

For the third example concerning the influence of atmo- 
spheric velocities on line profiles we refer to Sects. lBT2l and 
17.41 where we discuss the effects of a shock wave propagating 
through the line forming region of CO Av = 3 lines. This may 
not only result in asymmetric line shapes, but also in the phe- 
nomenon of line doubling for certain cases (Fig. 122b. 

For the fourth example we now turn to low-excitation CO 
Av = 2 lines and the example which can be found in Fig. 3 of 
Paper I. Compared to the second overtone CO lines, which orig- 
inate in a relatively narrow depth range, these first overtone CO 
lines are formed throughout the dust-forming region with con- 
tributions from many layers with diverse velocities (Paper I). 



0.8 1.0 1.2 1.4 

radius R [R. = 355R ] 

Fig. 11. Demonstration of the influence of atmospheric veloci- 
ties on line shapes. Upper panel: Gas velocities of one selected 
phase of the dust-free model W (black), overplotted is the cor- 
responding density structure (dash-dotted). Lower panel: The 
insert depicts synthetic CO Av = 3 line profiles calculated with 
(black) and without (grey) taking velocities into account for the 
radiative transfer. The main plot shows the radial optical depth 
distributions (central beam) for three points A/B/C in the line 
profile corresponding to different velocity shifts. The correlation 
is provided by the colour-code. 

Therefore the observed and the synthetic spectra show com- 
plex, multi-component line profiles, as discussed in detail in 
Papers I+II. Furthermore, the redistribution of opacity caused by 
atmospheric velocity fields may also influence the region of line 
formation of the strongest component. This was exemplified in 
Sect. 5.3 of Paperl for the results of the synthetic CO 2-0 R19 
line profile for one phase of model S. 

5. Simulating line profile variations 

It turned out in the past that high-resolution (near-IR) spectrosp- 
copy and the study of line profile variations represents the most 
promising tool to trace the kinematics in AGB atmosphere (e.g. 
GH04; their Sects. 4.8.3.1-2). This is exemplified in Fig.fHlwith 
the help of modelling results. 

5. 1 . Probing different atmospheric depths 

In their Sect. 4.8.1, GH04 pointed out that the atmosphere of 
an evolved AGB star (pulsating and mass-losing) can be sub- 
divided into three distinctive regions: the deep and dust-free 
photospheric layers dominated by the pulsating stellar interior 
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Fig. 12. Dynamics within the very extended AGB atmospheres 
results in complex, time-dependent structures (Sect. l2.2l ). This is 
illustrated in the upper panels, showing the corresponding gas 
densities, temperatures and velocities as well as the degree of 
dust condensation for one phase of model S. Together with the 
diverse properties of spectral lines (Sect.fJJ and the significant 
influence of velocity fields (Sect.|4]l, this leads to a complex line 
formation process. Thus, the synthetic line profiles of selected 
CO lines of the Av = 3 / 2 / 1 vibration-rotation bands computed 
for the respective model differ strongly, depending on the atmo- 
spheric regions where the lines are formed. This is demonstrated 
in the lower panels where typical lines (taken from Paper I and 
Fig. [131 are plotted for the cases that velocity effects are ac- 
counted for in the radiative transfer (black) or not (grey). The 
resulting profiles of different lines can be used to probe gas ve- 
locities at separated atmospheric layers (cf. Sect. 16. 21 Fig. 9 in 
Paper I). 



(zone A), the region of dust formation where the stellar wind is 
triggered (zone B), and finally the region of steady outflow at 
the inner circumstellar envelope (zone C). These differ strongly 
from the dynamic point of view. Our dynamical model atmo- 
spheres show a similar behaviour, as can be recognised in Fig.[T2l 
or Fig. |2] 



Observational spectroscopic studies (e.g. HHR82) demon- 
strated that individual spectral lines exhibit quite different line 
forming regions and can be used to gather information (e.g. ve- 
locity fields) concerning the above mentioned zones. Molecular 
lines of CO in the NIR proved to be especially useful in this 
context as discussed in Sect. 2.2 of Paper I or Sect. 1.3.2 of N05. 
A notable amount of high-resolution spectroscopy data of LPVs 
was collected during the past thirty years, basically by Hinkle, 
Lebzelter and collaborators. For a detailed review of the re- 
sults and a comprehensive compilation of references we refer 
to Sect. 1.3 of N05. From these investigations, it is known that 
second overtone (Av = 3), first overtone (Av = 2), and fundamen- 
tal (Av = 1) vibration-rotation CO lines roughly probe GH04's 
zones A, B, and C. Hinkle and collaborators derived excitation 
temperatures for certain types of lines by curve-of-growth anal- 
yses. The estimated temperatures in the regions of line forma- 
tion for CO Av = 3 / 2/ 1 lines amount to approximately 2000- 
4500/ 800-1200/ 300 K. 

In Papers I+II we studied the line formation including ve- 
locity effects for various molecular features (cf. Table|2| on the 
basis of the same dynamic model atmospheres as used in this 
work (Sect.[2]i. This led to the conclusion that the models from 
DMA3 are able to reproduce the finding that various spectral 
features originate at different geometrical depths within the at- 
mosphere. We estimated the temperatures in the line-forming 
region of CO Av = 3/2/1 lines to be about 2200-3500/800- 
1500/ 350-500 K (Paper I), which agrees well with the values 
derived from observations. Sampling atmospheric layers com- 
parable to the three zones specified by GH04, the synthetic line 
profiles of the different CO lines show quite different shapes if 
velocity effects are taken into account in the spectral synthesis. 
Figure[T2] illustrates the strong influence of atmospheric veloci- 
ties at different depths with characteristic CO line profiles from 
our model calculations. It also demonstrates how the gas veloc- 
ities convert into radial velocities, measured by the wavelength 
shifts of the deepest points of the complex line profiles. 

5.2. Tracing atmospheric movements over time 

Repeatedly studying one selected line profile in (observed) spec- 
tra and deriving RVs allows us to investigate the movement of 
the atmospheric layers where the respective feature is formed. 

A prominent example are the often used second overtone 
vibration-rotation CO lines, which can be found in NIR spectra 
at A « 1 .6 fim. As they are well suited to sample the inner photo- 
spheric layers where the movements are ruled by the pulsation, 
CO Av = 3 lines show a quite typical behaviour (HHR82) for 
Mira variables^ This is exemplified in the left panel of Fig. [13] 
with a time-series of line profiles (average of 10-20 unblended 
lines) of the S-type Mira^Cyg for various phases </> v during the 
lightcycle. Repeating in the same way every pulsation period, the 
following characteristic pattern is foundQ A blue-shifted com- 
ponent becomes visible before light maximum (0 v a:O.8). With 
increasing intensity it moves towards the red, crossing RV=0 
at around </> v ~0.4. Then the component becomes weaker again 

6 Note that there is a fundamental difference between Miras and 
SRVs, which is the topic of Sect.[7] 

7 The simplified picture of independently evolving components may 
be misleading to some extent, though, as the spectral features are 
formed over a whole region in depth and the line profiles (sum of dif- 
ferent contributions) can appear more complex (e.g. middle panel of 
Fig. l 13l >. Still, it can help to understand the scenario and for the certain 
phases of line-doubling the velocity jump due to the shock front results 
in two well-pronounced absorption components. 
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and disappears red-shifted shortly after the next light maximum 
(0 v ~O.l). This behaviour leads to the characteristic discontinu- 
ous, S-shaped RV curve as shown for example in Fig. [14] The 
occurrence of line doubling at phases around the visual maxi- 
mum is usually interpreted by shock waves passing through the 
line-forming region (Fig. 1 of Alvarez etal. BUOOl HHR82, N05). 
Based on a thorough analysis of these CO Av = 3 lines (RVs, ex- 
citation temperatures, column densities), HHR82 were able to 
derive a quite detailed picture of the photospheric structure being 
affected by the propagating shock wave as well as of the prop- 
erties of the shock front itself (emergence, progress, velocities, 
etc.) in their Sect. Via. 

In Fig.[13]we compare our modelling results (C-rich model) 
with observed spectra of the well studied S-type Mira ;f Cyg, 
as the spectroscopic monitoring of such CO Av = 3 lines for 
C-Miras is hampered by several problems (severe contamination 
by mainly CN and C2 in the H-band, availability of targets and 
instruments for long time series, etc.). Assuming that the princi- 
pal behaviour should be similar for all spectral types (because of 
the unique role CO; cf. Sect. 4.4.5 in GH04), such a comparison 
seems justifiable. In contrast to the results presented in Paper I, 
we neglected the contribution of other molecular opacity sources 
(C-bearing molecules, see Paper I for details) for the synthetic 
line profiles shown in Fig. [13] for an easier comparison. As de- 
scribed in detail in Papers I+II, we succeeded to reproduce the 
observed fact that CO Av = 3 lines trace the inner pulsating lay- 
ers of the atmosphere by our modelling approach. Nevertheless, 
there remained some shortcomings concerning the synthetic line 
profiles based on model S. For example, the components of dou- 
bled lines are not as clearly separated as in the observations. We 
also find quite complicated line shapes (at least for the highest 
spectral resolutions) during phases where observations show a 
transition of one component from blue- to red-shifts (e.g. phases 
V =O.57 and <p\, o \=0.33 in Fig. [Lit. This was the motivation for 
further efforts, the outcome of which will be discussed in Sect.[6] 

6. Towards realistic models for Miras 

In Papers I+II we were able to show that the used dynamic model 
atmospheres allow us to qualitatively reproduce the characteris- 
tic behaviour of spectral features sampling different atmospheric 
regions by consistent calculations. However, we also pointed out 
(Sect. 6 in PapeiTI) the shortcomings of the used model S and 
listed aims for subsequent modelling efforts (reasonable velocity 
amplitudes in the deep photosphere or global velocity variations 
in quantitative agreement with observed stars). Below, we try to 
vary the input parameters of the model (Table[TJ to get even more 
realistic atmospheric structures. The results of this preliminary 
parameter study is the topic of the next two subsections. 

6.1. The pulsating layers (CO Av = 3 and CN lines) 

The first and foremost aim for the fine-tuning of model parame- 
ters was a dynamic model that exhibits an amplitude of the gas 
velocities within the pulsating layers that is closer to what is 
found in observations. 

Hinkle et al. dl984t HSH84) state that though Miras are 
"clearly individuals" concerning the behaviour of CO Av = 3 
lines (differences in phases of line-doubling, range and phase- 
dependency of r exc , shapes and scales of column densities, dif- 
ferences in total RV amplitudes), they still share some charac- 
teristic features. The most important one is the behaviour of 
the line profiles as sketched in Sect. l5.2l and Fig. [13] resulting 
in the common shape of the velocity variations. Lebzelter & 
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Fig. 13. Left: Time series of average line profiles of CO Av = 3 
lines in FTS spectra of ^fCyg (cf. Fig. 1 in HHR82), taken 
from Lebzelter et al. (120011 ) who analysed the radial velocities. 
Observed heliocentric velocities were converted to systemic ve- 
locities RV* of the star assuming a CMRVheiiocent of -7.5 km-s -1 
(HHR82). Centre and right: Synthetic second overtone CO line 
profiles (5-2 P30) seen with spectral resolutions of 300 000 and 
70 000 for selected phases during one lightcycle, calculated on 
the basis of model S. 



Hinkle (2002a) presented a compilation of measured RVs for 
most of the Miras studied at that time (by Hinkle, Lebzelter 
and collaborators) in one plot, which has been adopted here 
in Fig. [14] It appears that Miras have a rather universal RV 
curve with discontinuities around maximum phases and a part 
where velocity increases linearly from blue- to red-shifts (neg- 
ative to positive RVs) through minimum light. The RV am- 
plitude (difference between minimum and maximum values) 
amounts to A/?Va:20-30km-s _I for all objects. This uniform pic- 
ture seems to be valid over a wide range of effective temper- 
atures, periods, (probably) metallicities and for different atmo- 
spheric chemistries. The last one was actually found only for 
M-/S-type stars, as measuring CO Av = 3 lines in carbon star 
spectra proves to be challenging (Barnbaum & Hinkle 1995). 

As it is a common feature of all Miras, realistic model cal- 
culations need to reproduce the discussed RV curve of Fig. [14] 
Our original model S was able to fulfill this qualitatively (cf. 
Fig.[T5ll. However, the step in gas velocity at the location of the 
shock was too small. Therefore, the splitting of CO Av = 3 lines 
appeared too weak (<pb o \=0. 75 in Fig.[T3l and the derived veloc- 
ity amplitude was too low (ARVa;14km-s _1 ). Model M, the re- 
sult of our small parameter study, was tuned to fit this aspect. 
Among the several modified parameters the increased piston ve- 
locity amplitude Am p is of particular importance. This led to a 
more realistic variability of the atmospheric structure in the in- 
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Observed radial velocities, Miras, CO Av=3 
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Fig. 14. Compilation of RVs derived from second overtone CO 
lines in FTS spectra of a large sample of Miras, demonstrat- 
ing the rather uniform behaviour. Observed heliocentric RVs are 
converted to systemic velocities by the respective CMRV of each 
object. The grey line was fit to guide the eye. Data adopted from 
Lebzelter & Hinkle ([2002a]). 



ner parts. The difference between post-shock outflow velocity 
and pre-shock infall velocity can reach values close to the esti- 
mateof »34km-s~ 1 given by Scholz & Wood (2000; summarised 
in Paper II), as can be seen in Fig.|4]for <f>b o i=0.72. The synthetic 
CO Av = 3 line profiles look similar to those shown in Fig.[T3l 
the RV curve - shown in Fig. [15]- appears much more like the 
observational results in Fig. [14] thoughQ Model M is able to re- 
produce this fundamental characteristic of Miras even quantita- 
tively. Although we needed to apply an arbitrary phase shift of 
Ac/>=0.3 to the synthetic RV curve in Fig. [15] to achieve agree- 
ment in phase with the observations in Fig. [14] the discontinu- 
ous RV curve resembles observations concerning the S-shape, 
the line doubling interval, the zero-crossing phase, the asymme- 
try w.r.t. RV = (infall velocities larger than outflow velocities as 
observed for most Miras), and - mainly - the velocity amplitude 
ARV of «21kms -1 . 

Hinkle & Barnbaum (119961 HB96) demonstrated how to cir- 
cumvent the difficulties in observing CO Av = 3 lines in C star 
spectra by investigating spectral features of CN in the NIR. Also 
sampling the inner photosphere of a Mira, they show a similar 
behaviour as second overtone CO lines (Sect. l5.2l Fig. [TBI. This 
can be recognised from the line profiles (left panel of Fig.lT6t. 
which repeat in the same way every lightcycle and the derive- 
able RV curves (right panel of Fig.fTTTi. On the basis of the im- 
proved model M, we synthesised CN line profiles following the 
approach of Paper! The results of the modelling for selected 
phases are shown in the middle and right panel of Fig. [16] It be- 
comes immediately apparent that compared to the profiles cal- 



8 Keeping in mind that the line profiles as well as the RVs may 
become more complex for the highest resolution of R = 300 000 (see 
Papers I+II), we will use only the RV results of the lower resolution of 
R = 70 000 to be compatible to observational results. 
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Fig. 15. Radial velocities derived from synthetic line profiles 
(spectral resolutions of R = 70 000) of second overtone CO lines 
calculated on the basis of model S (Fig. [13] cf. Figs. 2+9 in 
Paper II) as well as the corresponding ones for the new model M. 
An arbitrary shift of 0.3 in phase 0boi was applied for easier com- 
parison with the observational results in Fig. [14] which are given 
in c*v The typical S-shaped, discontinuous RV curve as well as 
the amplitude ARV found in observations (Fig.[T4l> can be repro- 
duced with model M. 



culated with the previous model S (see Fig. 5 in Paper I), the 
new line profiles reproduce observations much more realisti- 
cally. Although substructures may be identified at the highest 
resolution of /?=300000 (as for model S), the characteristic pat- 
tern of lines sampling the deep pulsating layers is clearly visi- 
ble, especially when rebinned down to the lower spectral reso- 
lution of 70 000. The transition of one component from blue- to 
red-shift (0boi~O.9-^-> 0.5) is much smoother. Pronounced line- 
doubling is found due to the more extreme velocity gradient 
across the shockfront. The splitting into two well separated com- 
ponents during certain phases can be recognised in Fig. [16] for 
<pbo\-0.72 (compare the corresponding velocity structure of the 
model in Fig.[4]i. The more realistic atmospheric structure of 
model M is also reflected in the RVs derived from the synthetic 
CN lines. The resulting discontinuous RV curve - included in 
Fig-US] - shows the same S-shape0 as the corresponding one 
from CO Av = 3 lines. The small phase shift between these two 
RV curves suggests that the line-forming region of CN lines lies 
somewhat further away from the centre of the star and the shock 
wave passes through at a later phase of the pulsation cycle Fl The 



9 Interesting is that this is not clearly recognisable in observed CN 
velocities of S Cep as shown in the upper right panel of the same figure. 
There the RV curve appears more like a straight line. One may suspect 
that the cross-correlation technique, applied by HB96 to obtain these 
velocities from the FTS spectra, is not able to resolve the weak com- 
ponents, as they appear shortly before light maximum and disappear 
shortly after the next maximum. 

10 The mean difference between the RV curves of the two types of 
lines (A0boi~O--/) together with an estimate of the propagation velocity 
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Fig. 16. Left: CN Av = -2 lines profiles (1-3 P 2 38.5) in the FTS 
spectra of S Cep observed by HB96 for several phases (cf. Fig 3 
in Paper II). Observed heliocentric velocities were converted to 
systemic velocities RV* of the star assuming a CMR Vheiiocent of 
-26.7 km-s - ' (HB96). Centre and right: Sequence of synthetic 
profiles of similar CN lines (1-3 Ch4.5) for different spectral 
resolutions of 300 000 and 70 000 throughout a lightcycle com- 
puted by using model M. 



velocity amplitude for the synthetic CN lines of ARV - ~22 km s 
resembles nicely the value of 22.3 found for S Cep by HB96 
(previously, we found «13.5kms _1 for model S). The above 
mentioned arbitrary phase shift of A0=O.3, which is needed to 
align the synthetic RV curve with the observed one in Fig.fTTl ap- 
pears plausible here as well. It is remarkable that the RV curve of 
S Cep is extending to more negative velocities w.r.t. the CMRV, 
while the modelling results are - in agreement with observa- 
tional as well as synthetic CO Av = 3 velocities - shifted to more 
positive values. 

Note that the stronger piston of the improved model M 
(A(/p=6kms~ 1 instead of 4kms _1 for model S; cf. TableQ} re- 
sults only in slightly higher outflow velocities, while the maxi- 
mum infall velocities are increased considerably, as can be seen 
in Fig. [15] The reason for this effect is the larger surface gravity 
(logg*) of model M compared to model S. Thus, the increase in 
total velocity amplitude ARV is mostly due to the faster infalling 
material. 

6.2. The global velocity field 

Simultaneous spectroscopic monitoring of various spectral fea- 
tures originating in different layers and for several instances of 
time (e.g. during a pulsation period) enables us to trace the evo- 



of the shock wave (;/f lont ~ 10. 16 km s '; see Sect.[8]( would result in a 
difference AR of approximately 0.15 R* = 61.8 R B = 0.288 AU. 



lution of the global velocity field throughout the outer layers of 
an AGB star and thereby the mass loss process. 

Extensive time series IR spectrocopy, which is needed to 
characterise the overall atmospheric dynamics, is available only 
for a few stars (limited availability of high-resolution IR spectro- 
graphs, demanding observations over years due to the long peri- 
ods). Investigated with sufficient coverage in phase and spectral 
range (i.e. line types) were mainly the O-rich Mira RLeo, the 
S-type Mira *Cyg, and the C-rich Mira SCep (N05). HHR82 
presented a fundamental study on x Cyg- Their time series of 
NIR spectra and RVs derived from different spectral features re- 
vealed remarkable details about the atmospheric motions of this 
starQ Although dealing with an S-type star, the results can be 
considered to be representative for Miras in general due to the 
singular role of the CO molecule (e.g. Sect. 4.4.5 in GH04) and 
the resulting presence in atmospheres of objects of all spectral 
types. Obtained radial velocities from HHR82 are compiled and 
plotted in the left panel of Fig. [17] supplemented with measure- 
ments of the expansion velocity of the circumstellar envelope 
(CSE) for comparison. There is only one C-type Mira, namely 
S Cep, for which a similar plot can be produced. The right panel 
of Fig. [17] shows the corresponding RVs compiled from HB96 
and Barnbaum d 1992bb . as well as again the range of CSE out- 
flow velocities coming from radio observations of molecular 
emission linesQ While CO Av = 3, CN and CO Av = 2 high- 
excitation lines sample the regularly pulsating layers of the deep 
photosphere, CO Av = 2 low-excitation lines are formed in the 
dust-forming layers. CO Av = 1 and selected atomic lines (e.g. 
KI at 7698.96 A) probe the outermost wind region. 

In Fig.[l8]the RV results of our line profile modelling based 
on model M are shown, aiming also for a global mapping of 
the dynamics throughout the atmosphere. As already described 
in detail in the previous section, the new model M proves to be 
rather realistic concerning the velocity variations in the deep at- 
mospheric layers (governed by the periodic pulsations) and the 
resulting line profile variations of molecular features originating 
there (i.e. CO Av = 3, CN). 

Observed low-excitation CO Av = 2 lines were found to 
probe the dust-forming region (Paper I, Sect. 15. 21) and show 
broadened, asymmetric line shapes (e.g. Fig. 3 in HB96). The 
multi-component profiles appear to be blends of a few contri- 
butions, which are at the highest resolutions currently avail- 
able not separable individually though (no pronounced line split- 
ting), and the temporal variations are not related to the lightcycle 
(cf. Papers I+II). Although all investigated Miras have in com- 
mon that the RV (measured from the deepest point of the main 
component) variations of these lines are not at all periodic (at 
least on the time scales of the observations), the behaviour can 
be somewhat different depending on the object considered. For 
X Cyg the RVs of CO Av = 2 low-excitation lines always stay 
close to the CMRV (Fig.fTTli. For other stars, as for example 
RVir (Lebzelter et al. 1999) or SCep (Fig.[T7), the velocities 
show moderate variations around the CMRV with amplitudes 
of up to a:15km-s _I (in any case lower than the amplitudes of 
CO Av = 3 lines though). Line profiles of the rather special star 
IRC +102 16 (dense and optically thick dust shell resulting from 
pronounced mass loss) are clearly blue-shifted indicating only 
outflowing material (Fig. 9 of Winters et al. 2000). These vary- 
ing behaviours may result from the fact that velocities in the 



11 ^fCyg appears to be the most extensively studied LPV from the 
kinematics point of view, at all, as outlined in Sect. 1.3.4 of N05. 

12 A similar compilation of RV data for the M-type Mira RLeo (col- 
lected by K. Hinkle and collaborators) can be found in Fig. 1.19 of N05. 
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Fig. 17. Radial velocities derived from (components of) observed or synthetic line profiles providing information on the gas veloci- 
ties in the respective line-forming region. 

Left: Radial velocities vs. visual phase during a lightcycle (<p v -Q corresponds to phases of maximum light in the visual) for several 
types of CO lines in observations of ^ Cyg. The data are adopted from HHR82. Observed heliocentric velocities are converted 
to systemic velocities of the star assuming a CMR Vheiiocent of -8.08 kms 1 , which results from applying the general conversion 
fiVheiiocent=^V L SR-17.98 (HHR82) to the CMRV LS r given by the more recent study of KYL98. Also marked (hatched blue) is the 
range of expansion velocities of the outflowing circumstellar material as measured by several publications on radio observations 
(DRM78, KYL98, NKL98, BGP89, LF093, RSOL06). 

Right: Radial velocities vs. visual phase for different lines in observations of S Cep. The data are adopted from HB96 and Barnbaum 
d 1992bb . Observed heliocentric velocities are converted to systemic velocities of the star assuming a CMR Vheiiocent of-26.65km-s~ 1 , 
which results from applying the general conversion /?Vheiiocent=-RVLSR-l 1-65 (HB96) to the CMRVisr given by BGP89. Again, the 
range of expansion velocities of the outflowing circumstellar material as measured by several publications on radio observations 
(B92a, NKL98, BGP89, LF093, OLN98, OEG93, RSOL06) is marked (hatched blue). 



line-forming layers are rather sensitive to the respective stellar 
parameters and the wind-acceleration process. In addition, opti- 
cal depth effects will have some influence. Synthetic CO Av = 2 
line profiles based on Model M show - apart from some tem- 
porally varying asymmetries, which can be interpreted as minor 
photospheric contributions - a strong main component, which 
appears blue-shifted at all phases. The slow but quite steady out- 
flow velocity of a:-5km-s _1 (Fig.[T8l) points towards a line for- 
mation in layers where the atmospheric material is already ac- 
celerated and variations due to shock waves are low (Fig.l4l)P^I 
In principle it could also be that the onset of the stellar wind in 
model M occurs in a too smooth way and the velocity variations 
in the dust-forming region are too small compared to a Mira like 
S Cep, which seems unlikely for model M, though (cf. Fig.@). 

However, general statements on the dynamics in the dust- 
forming regions of Miras based on high-resolution spectroscopy 
appear to be not within reach at the moment due to limitations 



13 A re-calculation of model M with an improved version of our RHD 
code - increased number of wavelength points for the RT, updated stick- 
ing coefficients, grain-size dependent treatment of dust instead of small 
particle approximation, etc. - led to a faster (u » 10-15 kms -1 ) but less 
dense (M~ 10~ 6 ) wind with a decreased dust content (f c xQ. 1). One 
may suspect that the decreased gas/dust absorption will lead to CO 
Av = 2 lines originating closer to the centre of the star, where the ve- 
locity variations are more pronounced. With the inner regions receiving 
the same mechanical input and a different wind structure at the same 
time, the line profiles based on this new model should therefore look 
more realistic. 



on the observational as well as on the modelling side. For the 
former we are confronted with a very small number of observed 
targets (cf. Sect. 1.3.2.2 in N05), sometimes sparsely sampled in 
pulsation phase, showing a variety of behaviours (e.g. Lebzelter 
et al. |1999l . As such low-excitation CO Av = 2 lines are formed 
over a wide radial range within the stellar atmosphere, we face 
an intricate line formation process (influence of velocity fields, 
different origins) and, thus, complex line profiles. These are very 
sensitive to the atmospheric structure (p-T-u) and the occuring 
amount of dust, which are both strongly depending on the prop- 
erties of the star (stellar parameters, pulsation characteristics). In 
addition, there may be (phase-dependent) deviations from spher- 
ical symmetry (e.g. Woitke 2006a, Freytag & Hofner 2008) and 
various layers of different (projected) velocities contributing to 
the emerging spectrum. In contrast to the general behaviour of 
CO Av = 3 lines in Mira spectra (Sect. 16. Il l, a cycle-to-cycle or 
object-to-object comparison of first overtone lines does not de- 
liver such a uniform pattern. The individuality of different AGB 
stars is even more prononunced in the remarkable velocity varia- 
tions derived from spectral features in the 4 fim range (Lebzelter 
et al. 1200 11 1. We found similar effects by our modelling, al- 
though we only used three models to compute synthetic first 
overtone CO lines. While models S and M exhibit a rather steady 
dust formation and CO Av = 2 velocities with almost no tempo- 
ral variations, things become quite different for the model pre- 
sented in Sect. 6.1 of Paper II (designated model S 1 in N05).This 
model shows pronounced cycle-to-cycle variations of dust for- 
mation resulting in outwards propagating dust shells (cf. Fig. 10 
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Fig. 18. Radial velocities derived from various synthetic line 
profiles (spectral resolutions of /?=70000) computed on the ba- 
sis of model M. An arbitrary shift of 0.3 in phase was applied 
for easier comparison with the observational results in Fig. [171 

in Paper II or Fig. 3.17 in N05), and thus remarkably changing 
velocity fields in the region of line formation for this type of CO 
lines. It will need further systematic studies - covering a large 
number of observational targets as well as models over a reason- 
able parameter range - to end up with a more general picture of 
the behaviour of CO Av = 2 low-excitation lines and its relation 
to the corresponding parameters of the stars or models PI 

Although observational studies of CO Av = 1 lines are diffi- 
cult due to severe telluric absorption, these lines still proved to be 
useful, as they originate in the cool wind region of evolved red 
giants. Only very few spectroscopic studies (e.g. Bernat 1981) 
dealt with these fundamental mode lines, all of them finding 
clearly blue-shifted spectral lines (also with multiple absorption 
components or PCygni-type profiles). The results of the small 
time series of spectra for^- Cyg obtained by HHR82 are included 
in Fig. [17] showing similar blue-shifts for every instance of time. 
We find the same behaviour for the synthetic line profiles based 
on the improved model M. The line shapes look very similar to 
those based on model S (Fig. 8 in PapeiT), although with some 
asymmetric distortion resulting in RVs (derived from the deepest 
point of the absorption), which are higher than the actual out- 
flow velocities of the model. The derived RVs of »;-10km-s _1 
(Fig.[T8Tl compare rather well to the values of x Cyg and other 
Miras, while the C-type Mira S Cep appears to have an outflow 
velocity higher by a factor of «2, placing this star on the upper 
end of the velocity distribution for optical C stars (e.g. Olofsson 
et al. [1993l t j2004l Ramstedt et al. l2ot)6l 

Our line profile modelling has now reached a state where we 
can conclude that the dynamic model atmospheres show global 

14 Such an investigation aiming at reasonably sampled RV diagrams 
as in Fig.[T7]or Fig.Q~8] appears rather challenging, though, because of 
the limited access to suitable spectrographs end the efforts needed on 
the modelling side. 



velocity structures that are in general quite close to real pulsat- 
ing and mass-losing AGB stars. In particular, the characteristic 
and univsersal behaviour of lines originating in the pulsating in- 
ner photosphere (CO Av = 3, CN) can now realistically be repro- 
duced by our models. The next step would be to fit in addition 
the velocities in regions where other types of lines originate (si- 
multaneously with one model for a given object). However, such 
a tuning of the models to achieve a specific fit for a certain star 
appears to be difficult. Apart from the problems on the observa- 
tional side (insufficient determination of stellar parameters like 
luminosity, mass or C/O; sparsely available RV data from high- 
resolution spectroscopy) and the questionable relation between 
dynamic models and observed targets (cf. Sect. 3 in Paper II), 
the modelling itself sets limits to the tuning efforts. The quite 
time-consuming process from a given set of model parameters 
(Table[TJ to an RV diagram (Fig.[T8l will impede the fitting pro- 
cedure. Still, the outcomes of line profile studies can provide 
important constraints (in addition to photometry, low-resolution 
spectra, interferometry, etc.) for a detailed study of selected AGB 
stars. 

7. CO Av = 3 line profiles of SRVs and Miras 

7.1. Observed line profile variations of SRVs 

The behaviour of second overtone CO lines in spectra of SRVs 
is somewhat different from those of Mira variables CSect. l5.2T >. 
Time-series spectroscopy of a few selected objects and results 
on radial velocities were published by Hinkle et al. d!997l l. 
Lebzelter ( 1999) and others (cf. Tables 1.4 and 1.5 in Nowotny 
120051 ). It appears that there is no RV curve uniformly shaped for 
all SRVs as it was found for Miras. The velocity variations with 
time can be smooth, continuous and regular to different degrees. 
In a simple approach one would assume that irregular RV varia- 
tions reflect irregularities in the pulsation of the interior (which 
are also supposed to be the reason for the irregular elements in 
the light curves of SRVs) in combination with optical depth ef- 
fects. The RVs are not strictly periodic, velocity values do not 
repeat from one period to the next. That only Doppler-shifted 
profiles but no line doubling can be found at any phase repre- 
sents the clearest distinction from Miras. Semiregular variables 
have significantly smaller RV amplitudes (difference between 
maximum and minimum RV) than Miras, most objects studied 
show only variations in the range of a few km-s . Only very 
few exceptions show amplitudes of 10-15 km-s -1 (among them 
WHya, cf. Sect.O- Hinkle et al. ( fl997l ) state that all LPVs 
have random velocity variations of the order of a few km-s -1 
(as possible explanations they list large-scale convective phe- 
nomena, atmospheric turbulence, non-radial pulsation, random 
changes in the periods of the stars, etc.), "masking the intrinsic 
RV variations" due to pulsation in which we are interested here. 
As they have the same order of magnitude as the low ampli- 
tudes of SRVs caused by pulsation, these random RV variations 
are more relevant for such stars compared to Miras. From curve 
of growth analyses Hinkle et al. (1997) deduced that CO Av=3 
lines exhibit roughly constant excitation temperatures over the 
lightcycle (typically 3300-3400 K) with variations of less than 
200 K. Conspicuous is - and this is another important difference 
to Miras - that the velocity distributions of the SRVs studied by 
Hinkle et al. £[997) and Lebzelter (119991 ) are clearly asymmet- 
ric w.r.t. the systemic velocity RV - 0, but in the opposite way. 
Negative offsets of several km-s -1 from the CMRV were found, 
the second-overtone lines appear blue-shifted for most or even 
all of the lightcycle. This means that only outflow is observed, 
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while no infalling material is seen in the spectra. The authors dis- 
cuss possible explanations for this puzzling result. Among sev- 
eral reasons (long secondary periods, etc.), effects of convection 
were proposed as a solution. From observations of our Sun we 
know that the interplay between the intensity of granules and ve- 
locity fields can result in an overall blue-shift. On the surface of 
late-type giants large convective cells are expected to occur (e.g. 
Freytag & H6fner2008). This could provide an explanation for 
velocity asymmetries (convective blue-shift). The question was 
later resumed by Lebzelter & Hinkle d2QQ2b] ), who did not ob- 
serve time series of individual stars, but rather chose a statistical 
approach by obtaining data at two epochs for a quite substantial 
sample of SRVs. The resulting velocities, together with some 
older measurements from time series of individual stars can be 
found in their Fig. 2. Again, the majority of the measurements 
are shifted w.r.t. the CMRV. Nevertheless, Lebzelter & Hinkle 
rule out that variable convective cells on the surface are the sole 
reason for variability in these stars (some measurements show 
red-shifted lines with RV > CMRV), indicating that we observe a 
mixture of pulsation, convection and long period variations still 
to be clarified. Also temperature effects may play a role (only 
layers with certain temperatures lead to spectral lines, but may 
be observable only at certain phases). 

Although the asymmetric velocity distribution is not yet fully 
understood and no typically shaped RV curve for all SRVs was 
found, Fig. 2 of Lebzelter & Hinkle (200 2b) suggests the fol- 
lowing general behaviour (as pointed out by the authors). RV 
variations have on average amplitudes of 3^4km-s _1 with most 
positive velocities around phases of light minimum (c* v ~0.5) 
and most negative values around light maximum {(/> v xQ.O). The 
whole distribution appears blue-shifted by «1 km-s -1 . The veloc- 
ity variations may be rather continuous, and line doubling does 
not occurF^I 

It is interesting to note that the latter does not necessarily 
mean that there are no pulsation-induced shock waves. For some 
SRVs, as e.g. W Cyg in Hinkle et al. (|T997), the hydrogen lines 
of the Balmer series show up in emission during some light cy- 
cles (which is a clear indication for shocks; e.g. Richter et al. 
I2001kj2003l and the summarised references therein), even though 
no doubled lines can be found. We will show by means of model 
calculations that this behaviour can be reproduced by one of our 
dynamic model atmospheres. 

7.2. The semiregular variable WHya 

One star, namely the M-type SRV W Hydrae, is singled out here, 
as it will play a role for the following comparison with model W. 
Radial velocity variations derived from CO Av = 3 lines were 
analysed by Hinkle et al. (1 19971 1 and Lebzelter et al. d2005al >. 
This object appears somewhat strange concerning its attributes. 
Although some irregularities in the shape of the lightcurve and 
the visual amplitude were found, it appears to be more like a 
Mira variable from its light variation (P = 361 d , AV = 3.9 m ). This 
is also supported by its location on sequence C (fundamental 
mode pulsators) in P-L-diagrams (Lebzelter et al. 2005a). Based 
on the limited sample of measurements, Hinkle et al. specu- 
lated that also the velocity variations could be interpreted by 
a Mira-like behaviour. This was ruled out later by Lebzelter et 
al. d2005al l using their extended series of spectroscopic obser- 
vations. From the resulting plot of RVs versus phase, shown in 
Fig. [19] WHya has clearly to be assigned to the group of (large 
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Fig. 19. Radial velocities derived from Doppler-shifted CO 
Av = 3 lines in FTS spectra of WHya (left panel of Fig.l20l. 
Observed heliocentric velocities are converted to systemic veloc- 
ities of the star assuming a CMR Vheiiocent of 38km-s~' (Hinkle et 
al. U997l l. Data adopted from Lebzelter et al. d2005a). The grey 
line was fit to guide the eye. 



amplitude) SRVs. It shows a rather periodic, continuous and al- 
most sinusoidal RV curve. No line doubling could be detected. 
Still, the velocity amplitude of a:15km-s _1 is clearly larger than 
the average for SRVs. WHya exhibits a well pronounced in- 
trinsic RV variation, perceivable even on top of the (above dis- 
cussed) random variations. It was suspected by Hinkle et al. 
(119971 ) that this star could be in a transition from the SRV stage to 
the Mira stage, an evolutionary scenario that was introduced with 
the detailed studies of pulsation of AGB stars (cf. Sect.[T]i. Later, 
Lebzelter et al. (2005a) concluded that a higher mass might be 
the reason for the remarkable velocity behaviour (no line dou- 
bling and small ARV compared to typical Miras also pulsating 
in the fundamental mode). 

7.3. Synthetic line profiles and radial velocities 

While studying synthetic high-resolution spectra on the basis of 
different dynamic model atmospheres, we found that model W 
(described in detail in Sect. 12. 21 ) is able to reproduce the be- 
haviour of CO Av = 3 lines in spectra of SRVs - particularly 
of W Hya. This model of a pulsating atmosphere was computed 
with a rather low piston amplitude of Am p =2km-s~ 1 . No dust- 
formation takes place in the outer layers, and thus no stellar 
wind can develop. Also, WHya shows only a very low mass 
loss rate of M ~ 2- 10~ 8 M yr -1 (Hinkle et al. [l997l . It shall how- 
ever explicitly be noted that the stellar parameters of (the C-rich) 
model W were not chosen to resemble (the O-rich star) W Hya. 
Still, a comparison of this one certain aspect - the RV varia- 
tions of CO lines - will be made, as the model reproduces the 
behaviour characteristic also for all other SRVs observed so far 
(Sect.|TB- 

The resulting synthetic line profiles for a representative num- 
ber of phases during the lightcycle are plotted in the two right 



20 



W. Nowotny et al.: Line formation in AGB atmospheres including velocity effects 
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R = 300.000 R = 70.000 
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Fig. 20. Left: Time series of average line profiles of CO Av = 3 
lines in FTS spectra of the large-amplitude semiregular variable 
WHya, taken from Lebzelter et al. (I2005at who analysed the 
radial velocities (cf. Fig.[T9b. Observed heliocentric velocities 
were converted to systemic velocities RV* of the star assuming 
a CMflVheiiocrat of +38km-s~ I (Hinkle et al. [19971 . Centre and 
right: Synthetic second overtone CO line profiles (5-2 P30) at 
the two different spectral resolutions indicated, calculated on the 
basis of model W for selected phases 0boi- 
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Fig. 21. Radial velocities throughout the lightcycle as derived 
from synthetic CO Av = 3 profiles based on model W (right panel 
ofFig.|20]l. 



panels of Fig.[20] Whenever the corresponding profiles of the 
observed FTS spectra were available, they are also plotted for 
comparison in the left panel. The figure demonstrates that the 
CO Av = 3 lines only show a single blue- or red-shifted compo- 
nent as well as rather asymmetric line profiles at highest reso- 
lution (R = 300 000) due to the non-monotonic velocity field of 
the atmosphere (the extreme case of phase (f>b o \=0.48 was inves- 
tigated in detail in Sect. 14. 21 and Fig.fTTTi. RVs derived from the 
deepest point of the line profiles are drawn in Fig. [21] The plot 
shows that our model W is able to reproduce the principal be- 
haviour of semiregular variables as sketched in Sect. UJH (maxi- 
mum RVs around light minimum, no line doubling). Concerning 
the asymmetry w.r.t. RV — as well as the velocity amplitude 
(ARV ~ 8 km-s" 1 for model W), the synthetic RV curve is more 
similar to the observed data of W Hya in Fig.[19]than to the gen- 
eral behaviour of SRVs shown in Fig. 2 of Lebzelter & Hinkle 
(2002b), though. Only a small difference in phase may be sus- 
pected, so that the bolometric phases of the model lag behind the 
visual phases of the observations by 0.1-0.2. 

7.4. Shock waves and the occurrence of line doubling 

As discussed in Sect. 17.11 observational studies of semiregu- 
lar variables - in contrast to Miras - did not find doubled CO 
Av = 3 lines for any phase during the light cycle. This behaviour 
is somewhat astonishing as one would generally expect shock 
waves in pulsating atmospheres, and also Balmer emission lines 
were observed in spectra of selected SRVs without CO line dou- 
bling. However, we found the same result for line profiles based 
on model W as discussed in the previous section. In Fig. [22] the 
puzzling fact shall be explored with the help of our modelling 
tools. For selected phases of model W (left panels) and model S 
(right panels), we present in this figure the atmospheric structure 
(gas velocities and the corresponding densities) together with the 
resulting CO line profiles and the radial optical depth at signifi- 
cant wavelength points within the synthetic spectra (A/B/C). 

For the model without mass loss in the left panels, a 
shock front can be seen in the velocity structure around 1 R+. 
Nevertheless, only a blue-shifted and slightly asymmetric com- 
ponent with RV- -2.65 km-s -1 (B) from behind the front shows 
up in the line profile. Practically no red-shifted component is 
formed. The very weak one, which can be vaguely discerned in 
the synthetic spectrum at RV~ lOkm-s -1 (C), would certainly 
not be recognisable in observed spectra. The reason for this is 
the sharp drop in density in front of the shock wave where the 
necessary gas velocities would be available. The optical depths 
at wavelength points redder than the rest wavelength of the re- 
spective line are not large enough so that a component would 
become visible in the spectrumM Although model W shows 
infall velocities of up to * 14km-s _I , these are not translated 
into the spectrum as the related outer layers do not contribute 
to the formation of the spectral line. As discussed in Sect. 3.3.1 
of N05, the gas densities of the model needed to form an indi- 
vidual component of the used CO Av = 3 line have to be higher 
than logpg as «— 11.5g-cm -3 , which is not met for the phase of 
model W displayed in Fig. [22] 

However, this requirement is fulfilled for the phase of 
model S shown in the right panels of the figure. Due to the for- 
mation of dust and the resulting wind, the atmospheric structure 



16 Note: the optical depths at every wavelength / RV point of the spec- 
trum increases dramatically at the inner edge of the dynamical model 
and further on within the hydrostatic extension (see Sect.O because of 
the increasing continuum opacities, though. 
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Fig. 22. Illustration of the fact that the occurrence of a shock wave within the atmosphere does not necessarily lead to line doubling 
in spectra containing CO Av = 3 lines. Demonstrated by representative phases of the dust-free model W (left) and the wind model S 
(right). In the upper panels the atmospheric structures are plotted. The lower panels show the radial optical depth distributions at 
certain wavelength / RV points as marked with A/B/C in the insets containing the synthetic line profiles (with and without taking 
velocity effects into account for the radiative transfer). The correlation is provided by the colour-code. See text in Sect. 17.41 for 
details. 



differs somewhat in the region around 1 which is relevant 
for the line formation of CO Av = 3 lines. The layers in front 
of the innermost shock front show smaller velocities compared 
to model W, but the densities are enhanced by roughly one order 
of magnitude. Therefore also the infalling material contributes to 
the spectrum formation. The non-zero velocities of the model re- 
sult in a redistribution of opacity and decreased optical depth for 
the rest wavelength (A). A distinct blue-shifted component with 
RV=-5km-s- 1 (B) arises from a narrow region (marked with a 
blue diamond) behind the shock front located at 0.9/?+, similar 
to the previous case of model W. Moreover, a clear increase of 
optical depth in the layers in front of the shock (marked with a 
red parallelogram) can be recognised in the lower right panel for 
the radial plot corresponding to wavelength point C if velocity 
effects are taken into account. This leads to an additional red- 
shifted component at RV- 7.9 km-s -1 . Altogether, two separated 
components can clearly be seen in the synthetic spectrum based 
on this phase of model S. 

The phenomenon of line doubling is therefore not a matter of 
the occurrence of shock waves but rather depends on the ques- 
tion if there is infalling material with high enough densities to 
produce optical depths which are necessary to give rise to a sec- 
ond component in the line profile. By a comparison of Figs.Q] 
[3] and [20] it can be understood that the complex temporal vari- 
ation of the atmospheric structure of model W (due to the pul- 
sation) does not lead to a second line component at any time 



during the light cycle. For <f>\, o \-0.99 there is only outflow, and 
a blue-shifted component arises. At the later phase of (f>b o \=0.16 
the density gradient becomes less steep, and the gas velocities 
slightly vary around Okm-s -1 . As they are formed over a rel- 
atively wide range, the spectral line becomes asymmetric, but 
does not show a pronounced shift. By the time of <p\, Q \-0.48, 
the density structure shows a significant steplike decrease due 
to the emerging Shockwave. Still, all of the infalling layers con- 
tribute to the resulting line profile (cf. discussion in Sect. l4.2l and 
Fig.fTTTi. which appears red-shifted. After this, the outflowing 
layers behind the shock become more and more important due 
to the ever increasing densities. At the same time, the infalling 
layers in front of the shock loose their influence because of the 
dramatic decrease in density. At no phase the model shows a 
clear velocity difference AM gas together with relevant densities in 
the outer layers. Thus, no line splitting can be found, as exem- 
plified in Fig.[22]for the phase (pb o \=0.68. 



8. Remarks on different velocities 

We note that the derived ARV of the two components of a dou- 
bled line profile does not provide information concerning the 
propagation velocity of the shock wave Mf lon t, but reflects the 
difference in velocity Aw gas of the outflowing gas behind the 
shock front and of the infalling matter in front of it. This will 
be illustrated on the basis of model M. With the help of the 
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plot of moving mass shells in Fig. [2] the shock propagation ve- 
locity can be estimated to d/?/d^boi ~ 1.5 R+/ cycle leading to 
"front ~ 10. 16 km s -1 . On the other hand, a typical velocity dif- 
ference of the model at the location of the shock amounts to 
AMg as » 29 km s _1 (cf. Fig.|4j, if we take the line doubling phase 
4>bo\-0-72 (Figs.[T6rr[T8li as an example. 

However, the most negative values for RV derived from e.g. 
CO second overtone lines (Fig.lLSb provide a hint on the shock 
wave velocity Mf ron t- At least a lower limit can be given. For ex- 
ample, we find RV «-6.25kms _I for the blue-shifted compo- 
nent of such a line at phase <pb o \=0.93 of model M0 Multiplying 
this with a reasonable conversion factor p (Sect. 14. Il l leads to gas 
velocities close to the maximum outflow of -9.6 km s 1 behind 
the shock front (compare the velocity structure for this phase in 
Fig. [4]: with the shock located at * 1.257?*), the latter being in- 
terrelated With Mfront. 

Another type of velocity, which can be studied with the 
help of our models, is the movement of the formed dust 
shells. Considering the one recognisable at zz3-5R* in Fig. 1 of 
Paper I we can estimate d/?/d0boi~ 1-296 or a shell velocity of 
«sheii~10.5kms _1 . 

9. Summary 

The outer layers of evolved red giants during the AGB phase dif- 
fer strongly from the atmospheres of most other types of stars. 
Due to pulsations of the stellar interior and the development of 
a stellar wind, AGB atmospheres are far from hydrostatic equi- 
librium and the atmospheric structures show strong local vari- 
ations (spatially and temporally). This leads to a complex line 
formation process. Many layers at different geometrical depths 
and with various conditions (T, p, u) contribute significantly to 
the finally observable spectra, which are dominated by numerous 
molecular absorption features. 

Going to high spectral resolutions for studies of individual 
spectral lines, the complicated velocity fields become especially 
important. Relative macroscopic motions with gas velocites of 
the order of lOkms 1 heavily influence line profiles and their 
evolution in time. Doppler-shifts of spectral features (or com- 
ponents of a line) represent an indicator for kinematics in the 
atmospheric regions where the lines are formed. Thus, spectro- 
scopic monitoring of molecular lines originating at different at- 
mospheric depths allows us to trace the overall dynamics and 
explore the process of mass loss in AGB stars. 

On the basis of a typical atmospheric structure and by as- 
suming chemical equilibrium as well as LTE, we studied the 
influence of the parameters (E exc and gf) of spectral features 
on the respective estimated line-forming region within the atmo- 
sphere (Sect.O. This is of particular importance for the different 
rotation-vibration bands of CO, which play an important role in 
the course of kinematics studies. 

We proceeded with the second aspect of this paper, namely 
the modelling of line profiles with the help of state-of-the-art 
dynamic model atmospheres. For this purpose, we used models 



17 Note that Wf lont becomes higher than the initial pulsation velocity 
resulting from the mechanical energy input due to the piston, which is 
in the case of model M not larger than 6kms~ 1 (Table[T}. While the 
sound waves propagate outwards, they steepen (as described by GH04, 
see their Fig. 4.15) and become strong radiating shock waves. 

18 This is the RV value derived from the corresponding line profile at 
the highest spectral resolution of R = 300 000, rebinning the spectrum 
to 70 000 influences the line shape and leads to a slightly different RV 
compared to those plotted in Fig.l 151 



which represent a numerical simulation of the most widely ac- 
cepted (at least for C-rich stars) scenario of pulsation-enhanced 
dust-driven winds. Atmospheric structures at certain instances of 
time were used as input for the spectral synthesis, for which we 
assumed conditions of LTE as well as equilibrium chemistry. A 
spherical RT code was applied that accounts for velocity effects. 

Spectroscopic observations of Mira variables revealed quite 
characteristic line profile variations over the light cycle for 
spectral features probing the inner dust-free atmospheric lay- 
ers where the atmospheric movements are ruled by the pulsating 
stellar interiors. We were able to tune the input parameters of a 
dynamic model (model M) so that we could quantitatively repro- 
duce the - apparently universal - discontinuous RV curve (line 
doubling at phases of light maximum, S-shape, amplitude, asym- 
metry) of CO Av = 3 (and CN) lines, which can be interpreted as 
a shock wave propagating through the region of line formation 
(Sect. l6.ll >. To match observed RV curves of these lines with syn- 
thetic ones based on model M we needed to introduce a phase 
shift of A<f> » 0.3 between bolometric phases 0boi and visual ones 
<p v so that the former would lag behind the latter (in agreement 
with photometric results). 

In observed spectra of SRVs the same spectral lines show 
a different behaviour. Although there seems to be no gen- 
eral RV curve for this type of stars (as for Miras), a RV 
distribution peaking around phases of light minimum and no 
line doubling is commonly found. Synthetic line profile vari- 
ations (Sect. l7.3l > computed on the basis of a pulsating model 
atmosphere (model W; no wind) are able to reproduce this 
and compare well to observational results of the semiregular 
variable WHya (showing an almost negligible mass loss of 
Ma!2TO _8 M yr _1 ). The same model was used to investigate 
the phenomenon of line doubling (Sect. l7.4l i. We argue that a 
shock wave propagating through the line forming region may 
not be sufficient for the occurence of doubled lines in every 
case. Large enough optical depths in the infalling layers may 
be necessary for a second, red-shifted component to be visible 
in the spectra. This could be an explanation for the fact that for 
some stars we find observational evidences for shock waves in 
their atmospheres (emission lines), but do not observe doubled 
CO Av=3 lines. 

The connection between gas velocities in the line-forming 
region and RVs derived from Doppler-shifted spectral lines was 
examined by means of artificial velocity fields (Sect. 14. it . On the 
basis of a hydrostatic initial model and under the assumption of a 
constant value for M gas at every depth point of the atmosphere we 
found conversion factors p of w 1 .2-1 .5. We also point out possi- 
ble difficulties for relating measured RVs with real gas velocities 
(e.g. inferring the velocity difference across a shock front from 
doubled lines) due to the complicated line formation process. 

In summary the models studied in this work represent the 
outer layers of pulsating AGB stars with or without mass loss 
quite realistically. At least for the C-rich case, the dynamic 
model atmospheres approach quantitative agreement with ob- 
servations. One important aspect in the comparison of mod- 
elling results with observational findings are line profile varia- 
tions in high-resolution spectra. Agreement in the velocity vari- 
ations of various molecular NIR features sampling different re- 
gions within AGB atmospheres can provide information on the 
interrelation of pulsation and mass loss or set constraints on the 
acting mass loss mechanism. In addition, this may be a criterion 
to constrain stellar parameters of a model. 

Although selected observed dynamic aspects (CO Av = 3 
lines in Mira spectra) are reproduced very well by the atmo- 
spheric models and also the global velocity field resembles what 
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can be derived from observations of typical long period vari- 
ables, a fine-tuning of the model parameters to fit certain objects 
appears to be more difficult due to the complex dependence on 
several parameters. 

In contrast to the encouraging results for the second over- 
tone lines of CO in both Miras and SRVs, the case of CO first 
overtone lines is a lot less clear-cut both from the observational 
and the theoretical point of view. These lines are presumably 
formed within the dust formation and wind acceleration region, 
and they display quite different types of behaviour in the few 
well-observed stars, with a considerable spread in variability. A 
similarly varied picture - probably an indication of the dynam- 
ical complexity of the corresponding layers - emerges from the 
models that we have analysed in detail here and in earlier papers. 
In particular, model M (which reproduces the second overtone 
lines nicely) shows basically no variations of the first overtone 
lines. We are, however, confident that further tuning of stellar pa- 
rameters can produce a model which shows the typical behavior 
observed for second overtone lines and, at the same time, more 
pronounced variations in the first overtone lines. 

A detailed comparison of phase-dependent line profiles 
derived from models with observations is presently limited 
by uncertainties in the relation between bolometric and visual 
lightcurves. Bolometric phases - which are quite directly cou- 
pled to global atmospheric dynamics (pulsation and resulting 
shocks) - are a natural choice for discussing model features. 
Observations, on the other hand, usually are given in terms of 
visual phases, as bolometric variations are not directly accessi- 
ble. In principle, visual phases can be derived for the models, 
but they may depend strongly on details of dust processes 
(grain opacities can dominate in the visual range for C-stars, 
cf. AppendixlAl. which may have only a weak connection with 
the dynamical phenomena we want to probe by studying the 
line profiles. A more comprehensive discussion of photometric 
model properties and their observed counterparts will be the 
topic of a separate paper. 
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Appendix A: Bolometric vs. visual light variations 

A. 1. Observational findings 

A direct comparison of observed and synthetic line profiles (or 
derived RVs) is affected because different phase informations are 
used in the two cases. On the one hand, observations are usually 
linked to visual phases (p v within the lightcycle (e.g. AAVSO 
measurements). On the other hand, for the model calculations - 
a priori - only bolometric phases from the luminosity lightcurve 
(e.g. Fig. [3]) are known, which is itself determined by the inner 
boundary (cf. Sect. l2.2l ). A relation between the two kinds of 
phase information would be desirable for comparing the mod- 
elling results with observational findings. 

Measuring bolometric lightcurves of LPVs is not straigth- 
forward. Although a few authors estimated such variations from 
multi-epoch observations in several filters by fitting blackbody 
curves (e.g. Whitelock et al. 2000), a bolometric lightcurve is 
in general not available for a given star. However, it is known 
that lightcurves in the NIR (e.g. the photometric K-band) trace 
bolometric variations more closely than lightcurves in the op- 
tical wavelength range. The reasons for this are that (i) most 
of the flux is emitted in the NIR, and (ii) the visual spec- 
tra are severely influenced by features of temperature sensitive 
molecules Infrared lightcurves may therefore provide clues on 
the desired 4>bo\- <t>v - relation. 

Several studies compared variations in the visual and NIR 
of Mira variables in the past (cf. the list of references given by 
Smith et al. l2006l) . leading to a rather general result: the maxima 
of IR (i.e. ~bolometric) lightcurves lag behind the maxima of 
visual lightcurves by ^0. 1-0.2 in phase. A good illustration of 
this is provided e.g. by Fig. 2.44 of Lattanzio & Wood (2004), 
who show lightcurves of the M-type Mira RRSco in various 
photometric bands. By using a well chosen narrow-band filter at 
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Fig. A.l. Low-resolution opacity sampling spectra based on the 
initial model as well as selected phases of model W (dust-free, 
no wind). 



1.04yimi (located in a wavelength region with almost no absorp- 
tion features), Lockwood & Wing (I197U and Lockwood (I1972t 
also found a typical shift of A<f>x0. 1 if the phase values measured 
by this filter are assumed to correspond to bolometric ones. 
Recently, Smith et al. (2006 ) compared IR photometric data 



19 The latter effect is especially pronounced in the O-rich case and 
the characteristic bands of the fragile TiO molecule. It may be of less 
importance in the C-rich case as the species responsible for features in 
the visual, CN and C2, are more stable. 



( 1 .25-25ytzm) of the DIRBE instrument onboard the COBE satel- 
lite with visual observations (mainly AAVSO) for a sample of 
16 Miras and 5 SRVs. Their results confirmed the phase lag of 
0.1-0.2 between visual and IR phases for Miras. In contrast, 
they found no phase lags for SRVs. Smith et al. could repro- 
duce both results also with synthetic spectra based on dynami- 
cal model atmospheres of the Australia-Heidelberg collaboration 
(O-rich, dust-free). They conclude that phase lags are in general 
more likely or larger for more evolved stars (higher L, lower T e s, 
higher M, higher Amt, i, larger periods), that is for Miras rather 
than for SRVs. 

Most of the investigations in this context were carried out 
for M-type stars, results for C-rich LPVs are rather scarce in 
the literature. Thus, only some indications can be listed here. 
Kerschbaum et al. (2001 ; their Fig. 5) found that their lightcurve 
in the K-filter of the C-rich Mira TDra lags behind the visual 
data (AAVSO) by 0.06 in phase. Smith et al. (120061) state that 
they find no differences in phase lag for the spectral types M/S/C. 
For the one C-type Mira in their sample, V CrB, Smith et al. de- 
rive shifts of A0=O.13-O.15, the C -rich SRV UXDra does not 
show any shift. Smith et al. (2006) discussed optical-IR-offsets 
of broadband colours computed on the basis of the same type 
of C-rich dynamic model atmospheres as used here. They report 
that some of the models show phase lags, but do not give num- 
bers. However, no clear trend in the sign of the phase shift was 
found. For some models, visual phases precede IR ones, while it 
is the other way round for other models. 

A.2. Modelling results 

Following the suggestions of the referee, we studied the pho- 
tometric variations in the visual of the models used in this pa- 
per. Based on the dynamic model atmospheres of Tab.Q] and 
consistent with the descriptions in Sect. 12.41 we calculated low- 
resolution opacity sampling spectra (cf. Gautschy-Loidl et al. 
2004) with extensive wavelength coverage. By convolving these 
with transmission curves (Bessell & Brett 1988, Bessell 1990) 
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Fig. A.2. Same as Fig. lA. ll but for model M (dusty wind). 



for filters of the standard Cousins-Glass-Johnson system and ap- 
plying adequate zeropoints we get synthetic photometry. For de- 
tails of the method we refer to Aringer et al. (2009il. 

Example spectra (rebinned to R — 360) are shown in 
Figs. lA. lklA.2l demonstrating the qualitative difference between 
a model for a moderately pulsating red giant on the one hand 
and an evolved Mira with dust formation and mass loss on the 
other hand. In the first case, the spectrum shows only minor 
variations around the hydrostatic case during the pulsation cy- 
cle. The significantly changed radial structure in the latter case 
(Sect. l2.2l ) is also reflected in the spectra based on wind mod- 
els. Due to the characteristic absorption behaviour of dust grains 
composed of amorphous carbon (e.g. Fig. 1 of Andersen et al. 
1999), the whole spectral energy distribution is shifted towards 
longer wavelengths. It does not resemble the hydrostatic spec- 
trum at any point in time. 

For an easier comparison the synthetic lightcurves in vari- 
ous filters are plotted with fluxes normalised to [0,1] below|3 As 
shown in Fig. lA. 31 the light variation in the V-filter of the moder- 
ately pulsating model W follows the bolometric lightcurve quite 
closely, although some slight deviations from the pure sinusoidal 
shape are apparent. 

Figure lA~4l shows the light variations in various filters for 
the model M with mass loss. As found in observational stud- 
ies, IR-lightcurves stay rather close to the bolometric one. The 
lightcurves in the visual and red show some deviation, though. A 
sinusoidal behaviour with a phase lag of < 0. 1 is found for most 
phases throughout the lightcycle (ascending branch), with visual 
lightcurves lagging behind the IR/bolometric ones, in contrast to 
observational results. Significant differences between filters are 
seen on the descending branch of the lightcurve. The more or 
less pronounced asymmetries - monotonically increasing with 
decreasing wavelengths - suggest a relation with the absorption 
by amorphous carbon dust. 

Figure lA31 confirms this idea. If the dust absorption is not 
taken into account for computing spectra and photometry (not 
consistent with the dynamical modelling but a useful test), no 
asymmetry and phase shift is found for the V-lightcurve. If the 
dust absorption is included, the V-band flux of the star becomes 
considerably fainter and the asymmetry of the visual lightcurve 
becomes apparent. The pronounced decrease in V as seen in 
Fig. IA.5l is related to the new dust shell emerging just before light 
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Fig. A.3. Bolometric variation of model W (caused by the vari- 
able inner boundary; Sect. 12. 2t during one pulsation cycle and 
the resulting synthetic lightcurve in the visual. 
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20 Note that the amplitude of the variations may differ strongly be- 
tween filters, with amplitudes decreasing with increasing wavelengths. 
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Fig. A.4. Synthetic normalised lightcurves for model M for dif- 
ferent filters throughout one pulsation cycle. 



minimum for model M (see Fig.|2|. The increasing dust absorp- 
tion causes deviations from a sinusoidal variation and shifts the 
time of minimum (and maximum) light towards higher values 
of c*boi- This may provide an explanation for why a phase dis- 
crepancy of m 0.3 is found in Sect. l6.1l (Mira observations vs. the 
dusty model M), while only a value of as 0.1 is needed in Sect. 17.31 
(SRV observations vs. the dust-free model W). Assigning one 
global value for Acp = (pboi - <Pv is impeded by these deviations 
from a sinusoidal behaviour. 

Considering the possibly dominant role of dust opacities 
at visual wavelengths, however, any conclusions drawn from 
comparing observed and synthetic light curves (phases) should 
be regarded with caution for the current models. Recent mod- 
elling results indicate that grains in AGB star winds may grow 
to sizes where the small particle limit (used for the dynami- 
cal calculations here) is no longer a good approximation for 
dust opacities, in particular at wavelengths shorter than about 
1 micron (cf. Mattsson et al. 2009, Mattsson & Hofner in prep.). 



26 



W. Nowotny et al.: Line formation in AGB atmospheres including velocity effects 




phase c^ bol 

Fig.A.5. Synthetic lightcurves (absolute magnitudes) based on 
model M in the V-filter, where the absorption of circumstellar 
carbon dust is taken into account (same as Fig. IA.4b or not for 
the spectral synthesis. 

Unfortunately, model M is no exception in that respect, with a 
mean grain size of about 10~ 5 cm. For these grains, both the ab- 
solute and relative values of opacities at different wavelengths 
become dependent on the size of the particles. This may have a 
significant impact on photometric light curves. Phase-dependent 
effects due to grain growth may occur in different filters, pos- 
sibly affecting the relative positions of maxima and minima. A 
detailed investigation of this effect is, however, beyond the scope 
of this paper and will be discussed in a future paper about syn- 
thetic photometry. 

In summary, at the current state of modelling, bolometric 
light curves can be regarded as a quite reliable indicator of the at- 
mospheric dynamics of the models, since they are closely linked 
to the effects of pulsation and atmospheric shocks. In contrast, 
the diagnostic value of the synthetic V-lightcurves - and, con- 
sequently, of the synthetic visual phases - is rather doubtful for 
this purpose as the fluxes in the visual are strongly dependent on 
grain opacities (which probably require improvements in the fu- 
ture) and time-dependent dust processes, which may only have 
a weak connection with atmospheric dynamics. Therefore, with 
the additional uncertainties concerning the 0boi-0v- relation on 
the observational side (as outlined in the previous subsection), 
we decided to use the more representative bolometric phases 
when discussing modelling results in this paper. 



